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ABSTRACT 


The  P-3  Orion  patrol  airplane  has  a  netd  for  an  air-to-air  missile  system  for  defense 
against  enemy  aircraft  on  its  long-range  missions.  In  response  to  this  need,  the  Naval  Air 
Test  Center  was  tasked  in  1989  to  conduct  a  P-3/A1M-9  (Sidewinder)  integration  program. 
In  support  of  this  program,  a  \abration  test  stand  was  established  at  NPS,  and  a  ground 
vibration  characterization  was  conducted  to  determine  if  a  potential  flutter  problem  existed. 
This  test  resulted  in  the  development  of  a  twc  degree-of-freedom  lumped-mass  model  and 
experimental  determination  of  the  missile's  lesonance  modes  in  pitch,  With  the  ^ecen^ 
termination  of  the  P-7A,  the  P-3  community  is  now  looking  at  the  Orion  II  program  to 
carr\'  it  into  the  21st  century.  The  Orion  II  will  most  likely  have  a  beefed-up  wing 
structure,  necessitating  an  analysis  of  this  wing  in  conjunction  with  the  AIM-9  missOe. 
This  investigation  responds  to  that  requirement  by  concurrently  developing  a  mathematical 
model  of  the  AIM-9  missile  using  finite  element  methods  to  analytically  determine  its 
modal  parameters,  and  setting  up  a  modal  test  system  to  quantify  the  parameters  of  this 
model  by  experimentally  determining  the  inissile's  natural  frequencies,  mode  shapes  and 
transient  response.  This  fully  instrumented  test  system  and  associated  methodologies 
could  then  be  the  basis  for  conduf  ting  a  comprehensive  modal  test  of  the  AlM-9  missile 
system,  as  well  as  to  quantify  the  vibration  characteristics  of  other  candidate  missile 
sy.stems  for  the  P-3  and  its  eventual  successor. 
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INTRODUCTION 


1  . 

A .  BACKGROUND 

Over  the  last  25-30  years,  the  L'-3  Orion  patrol  airplane 
has  pel  formed  its  primary  mission  of  Anti-Submarine  Warfare 
(ASW)  with  virtually  no  defensive  weaponry  to  protect  itself 
on  its  solitary,  long-range  missions  on  the  world's  oceans. 
During  its  first  15-20  years  of  existence,  the  P-3  only 
carried  armament  in  the  form  of  mines  and  torpedoes,  as  well 
as  "Zuni"  rockets  for  gi ound  suppression.  It  had  been  viewed 
primarily  as  a  patrol  airplane  with  minimal,  if  any,  offensive 
capability.  However,  in  the  late  1970' s  and  early  1980' s, 
when  the  Harpoon  anti-ship  missile  was  integrated  into  the 
platform,  the  r-3  suddenly  became  an  offensive  weapon 
platform.  Throughout  all  this  time,  it  has  remained  highly 
vulnerable  to,  and  defenseless  against,  any  and  all  air 
threats . 

The  issue  then  becomes  one  of  increasing  the  survivability 
of  the  F-3  against  airborne  threats  .  During  the  last  few 
years.  Defensive  Air  Combat  Maneuvering  (DACM)  has  beeii,  and 
is  currently  being  taught,  to  train  P-3  aircrews  in  tactics 
that  would  lengtlien  the  P-3's  survival  time  in  an  air-to-air- 
engagement  .  This  training  has  shown  promise,  but  it  was  soon 
realized  that,  if  an  air-to-air  defensive  missile  were  to  be 
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iiiLoyi-aLed  iiiLo  Cho  weapon  invonl.ojL y ,  t-he  r-3  would  have  a 
nnich  better  chance  ot  set. vivinq,  it  not  possibly  destioying, 
the  ail  threat. 

This  realir.ation  led  1 1.'  tlie  r-3/AlM-9  pi ooi -ol -concept 
integiation  piograiu,  begun  in  1909,  and  currently  on-going  at 
t.lie  Naval  Aii  Test  Center  (NATC)  at  Naval  Air  Station  Patuxent 
Ivivei  ,  Maryland.  In  this  elfort,  NATC  was  tasked  to  evaluate 
the  separation  characteristics  o£  t.lie  AIM— 9  Sidewinder 
missile,  shown  in  Figures  1  ai'.d  2,  as  installed  on  P  — 3  outer 
wing  station  nine,  and  to  define  system  requirements  for 
pos.niblf'  use  in  a  P-3.  Wing  station  nine  is  of  critical 
iiuportanco  to  the  test  }irogram  because  it  allows  the  carriage 
of  an  AlM-9  missile  without  losing  a  Harpoon  wing  station. 


Figure  1  Sidewinder  Missile 

One  of  the  main  concorn.s  with  carrying  AIM-9  missiles  v'as 
whether  or  not  the  natural  vibration  modes  and  frequencies  of 
the  missile  sy’Stem  would  interfere  with  the  t:wo  low  trequeiK:y, 
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outer  wing  vibration  modes  of  the  P-3:  the  outer  wing  bending 
mode  at  4.7-8  hertz  (hz)  and  the  outer  wing  torsion  mode  at 
17-22  hz .  This  could  be  determined  by  conducting  a  Ground 
Vibration  Test  (GVT; ,  the  results  of  which  could  then  be 
incorporated  into  a  computer  model  of  the  P-3  wing,  leading  to 
a  determination  of  whether  or  not  a  potential  flutter  problem 
exists  with  an  AIM-9  installation  on  a  P-3  wing. 

In  support  of  this  on-going  program,  Hollyer  [Ref.  1] 
designed  and  constructed  a  GVT  stand  at  the  Naval  Postgraduate 
School.  A  ground  vibration  characterization  of  the  AIM-9 
missile  was  conducted,  resulting  in  the  development  of  a  two 
degree-of-f reedom  lumped-mass  model,  and  experimental 
determination  of  the  missile's  resonance  modes  in  pitch. 

With  the  recent  termination  of  the  P-7A  and  P-3H  programs, 
the  P-3  community  is  now  looking  at  the  P-3  Orion  II  program 
to  carry  it  into  the  21*’^  century.  The  Orion  II  will  most 
likely  have  a  beefed-up  wing  structure,  necessitating  an 
analysis  of  this  wing  in  conjunction  with  the  AIM~9  missile. 
This  investigation  responds  to  that  requirement  by 
concurrently  developing  a  mathematical  model  of  the  AIM-9 
missile  using  finite  element  methods  to  analytically  determine 
its  modal  parameters,  and  setting  up  a  modal  rest  system  to 
verify  the  mathematical  model  by  experimentally  determining 
the  missile's  natural  frequencies,  mode  shapes  and  transient 
response.  The  analytical  modeling,  combined  with  the  fully 


instrumented  experimental  modal  testing  pursued  in  this 
investigation  will  lead  to  the  following: 

•  Full  characterization  cf  the  dynamic  characteristics  of 
the  AIM-9  missile.  When  combined  with  an  appropriate 
analytical  model  of  the  wing,  these  methodologies  will 
allow  the  analytical  assessment  of  different  missile 
mounting  locations  and  configurations  without  the  risk  and 
expense  of  flight  tests. 

•  Modification  of  the  parameters  in  the  analytical  AIM-9 

model  to  analytically  assess  the  changes  in  the  modal 
characteristics  for  the  upgrades  and  modifications  of  the 
missile.  Upgrades  and  modifications  may  include 

propellent  change  in  the  rocket  motor  for  range  extension, 
and  modified  guidance  and  control  system  resulting  in  a 
mass  change.  Again,  this  is  possible  without  the  expenses 
of  prototype  construction  and  testing. 

•  Characterization  of  other  candidate  missile  systems  (e.g., 
AGM-65  Maverick  and  AGM-88  HARM)  for  the  P-3  and  its 
eventual  successor. 


B.  MISSILE  MODELING  TECHNIQUES 

To  obtain  a  structure's  vibration  characteristics,  two 
different  methods  can  be  used.  The  first  method  is  analytical 
modeling  which  utilizes  the  finite  element  method,  implemented 
in  a  computer  software  program.  This  is  a  "forward-type"  of 
analysis,  whereby  the  structure  is  partitioned  into  components 
whose  mass,  stiffness  and  damping  are  parameterized,  and  the 
structure' s  natural  frequencies  and  mode  shapes  are  predicted 
from  these  input  parameters.  The  second  method,  experimental 
modeling,  is  usually  referred  to  as  modal  analysis.  This  is 
an  "inverse-type"  of  analysis,  whereby  the  actual  structure  is 
subjected  to  a  known  forcing  function,  and  the  structure's 
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response  to  the  excitation  is  measured  at  multiple  locations . 
From  the  individual  responses  at  different  locations,  the 
overall  structure's  stiffness  and  damping  can  be  determined, 
along  with  the  natura],  frequencies  and  modes  of  vibration.  A 
block  diagram  of  both  of  these  methodologies  is  shown  in 
Figure  3  [Ref.  2:p.l-7]. 


Figure  3  Modeling  Techniques 


C.  IMPORTANT  DEFINITIONS 

A  compilation  of  terminologies  is  presented  to  facilitate 
the  discussions  of  the  finite  element  and  modal  analysis 
techniques.  Definitions  of  the  terms  are  available  from 
textbooks  of  vibration  analysis,  and  the  following  listing  was 
adopted  from  two  of  these  books  [P.efs.  3  and  4]  . 

• 


Modal  parameters  -  natural  frequency,  damping  and  mode 


•  Natural  frequency  -  a  frequency  at  which  a  system  vibrates 
when  disturbed  from  rest  and  then  released;  that  frequency 
where,  in  theory,  a  minute  force  can  produce  large 
movements  of  the  structure 

•  Resonance  -  occurs  when  a  force  coincides  with  the  natural 
frequency;  movement  of  the  structure  increases  either 
until  it  is  limited  by  damping  or  non-linear  effects,  or 
until  something  breaks 

f  •  Mode  -  resonance  or  peak  in  a  Frequency  Response  I’unction 

•  Modal  model  -  set  of  modal  parameters  representing  the 

*  dynamic  behavior  of  a  structure 

•  Modal  mass  -  mathematical  concept  that  can  be  thought  of 
as  being  proportional  to  the  amount  of  mass  actually  in 
motion  in  a  given  mode  shape;  it  is  not  the  physical  mass 
of  the  structure  and  cannot  be  measured 

•  Modal  stiffness  -  mathematical  concept  that  can  be  thought 
of  as  being  proportional  to  the  amount  of  stored  energy  in 
a  structure  in  a  given  mode  shape  at  the  instant  it 
reaches  the  point  of  extreme  motion 

•  Damping  -  quantity  that  controls  the  sharpness  of  a 
resonance;  usually  expressed  in  percent  of  critical 

•  Frequency  Response  Function  (FRF)  -  ratio  of  the  response 

(either  displacement,  velocity  or  acceleration)  at 
coordinate  i  to  the  excitation  at  coordinate  j  (for 
example,  =  X^/Fj) 

•  Measurement  -  FRF  acquired  frora  a  spectrum  analyzer 

•  Degree-of-f reedom  (DOF)  -  measurement  point  and  direction 
defined  on  a  structure 

•  Driving  point  FRF  -  measurement  made  with  the  excitation 
and  response  at  the  same  point  and  direction 

•  Reference  transducer  —  transducer  connected  to  the 

^  reference  channel  (channel  A)  of  the  analyzer;  normally, 

this  is  the  force  transducer 

t  •  Response  transducer  -  transducer  connected  to  a  response 

channel  (channel  B,  C,  or  D)  of  the  analyzer;  normally, 
this  is  an  accelerometer 
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•  Mode  shape  -  a  set  of  shape  coefficients  at  different 
locations  along  a  structure  that  defines  the  structure's 
deformation  at  a  specific  frequency 

•  Real  (normal)  shape  —  purely  real  (or  imaginary)  numbers 
that  describe  the  vibration  pattern  of  the  structure;  all 
parts  of  the  structure  are  moving  either  in  phase,  or  180® 
out  of  phase,  with  each  other 

•  Complex  shape  —  when  each  deformation  value  requires  two 
numbers  (magnitude  and  phase,  or  real  and  imaginary)  to 
describe  the  vibration  at  a  location;  can  have  any  phase 
relationship  between  different  parts  of  the  structure 

•  Curvefitting  -  matching  a  mathematical  equation  to  a  set 
of  data  points  obtained  from  an  FRF  measurement 
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IX.  ANALYTICAL  MODELING 


A.  FINITE  ELEMENT  METHOD  BACKGROUND 

In  designing  structures,  the  finite  element  method  is 
often  used  to  determine  the  structure's  natural  frequencies 
and  mode  shapes  before  it  is  actually  built,  in  an  effort  to 
foresee  any  potential  implementation  problems  which  may  occur. 
The  finite  element  method  can  be  viewed  as  a  building  block 
approach  whereby  the  stiffness,  damping,  mass  and  material 
properties  for  the  individual  building  blocks  are  known,  or  at 
least  estimated,  in  advance  and  coded  into  one  of  the  finite 
element  analysis  software  programs  currently  on  the  market. 
In  the  initial  design,  the  designer  will  have  at  least  a 
general  idea  of  the  configuration  of  the  structure,  the 
location  of  components,  and  any  mechanical  interfaces  and 
attachments.  Constituent  properties  of  the  structure  can  be 
estimated,  with  iteration  and  confirmation  of  the  changes 
leadj.ng  to  the  final  design  of  the  structure. 

The  finite  element  method  is  not  just  a  usefjl  tool  for 
designers .  It  can  also  be  used  to  determine  the  natural 
frequencies  and  mode  shapes  of  structures  which  are  already 
built.  Generally  speaking,  in  the  finite  element  method, 
structures  are  represented  by  discrete  grid  or  node  points 
which  are  connected  by  structural  elements.  In  other  words. 
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tlie  finite  element  method  can  be  viewed  as  a  discrete 


representation  of  a  continuous  system  made  in  order  to 
simulate  its  structural  behavior  and  response  to  expected 
forces  or  loadings .  This  representation  is  in  the  form  of  a 
mathematical  model  consisting  of  discrete  elements  connecting 
discrete  nodal  points.  The  more  node  points  that  are  chosen 
to  represent  the  structure,  the  more  accurate  the  solution 
will  be.  Conversely,  the  fewer  the  number  of  node  points 
chosen,  the  stiffen  the  model  will  be.  If  the  aim  of  the 
designer  or  analyst  is  a  detailed  stress  analysis  of  the 
structure,  then  the  nodal  density  must  be  increased  in  regions 
of  large  stress  gradients,  such  as  around  attachment  points 
and  around  the  location  of  applied  forces.  However,  if  the 
aim  is  a  deflection  shape  analysis  only,  then  fewer  nodes  may 
be  used.  [Ref,  5:Section  l,pp.  2-2,2-'3] 

For  either  an  existing  structure  or  one  that  has  yet  to  be 
built,  a  typical  finite  element  analysis  requires  the 
following  [Ref.  5;Section  l,p.  2-1] : 

•  Nodal  point  spatial  locations 

•  Structural  eleiaents  that  connect  the  nodal  points, 
representing  the  stiffness  of  the  structure 

•  Mass  properties  of  the  structure 

•  Boundary  conditions  or  structural  constraints 

•  Static  and  dynamic  load  application 
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The  mathematical  model  that  results  from  the  analysis  can  then 
be  used  to  determine  the  structure's  natural  frequencies,  mode 
shapes,  stress  levels  and  displacements.  An  added  benefit  of 
conducting  a  finite  element  analysis  on  a  particular  structure 
is  that  design  iterations  can  easily  be  accommodated,  cheaper 
and  faster  than  by  altering  the  actual  hardware.  Fo  example, 
if  the  designer  or  analyst  wished  to  determine  the  effect  of 
additional  mass  on  the  structure's  natural  frequencies,  all 
that  would  be  required  is  modifying  a  few  linos  of  software 
code,  and  rerunning  the  analysis.  The  effect  of  changing  the 
material  type,  damping,  shape  or  size  of  the  structure  can  be 
handled  easily  as  well.  However,  to  achieve  an  accurate 
model,  care  must  be  taken  in  the  assumptions  made  about  the 
structure's  geometry,  material  properties  and  so  on. 

B.  MSC/PAL2  FINITE  ELEMENT  ANALYSIS  SOFTWARE 
1 .  Introduction 

The  finite  element  analysis  software  package  used  in 
this  investigation  was  the  MacNeal-Schwendler  Corporation's 
MSC/PAIj2  Advanced  Stress  and  Vibration  Analysis  Software, 
Version  4.0.  MSC/PAL2  uses  the  finite  element  method  to  solve 
for  a  structure's  displacement,  forces  and/or  stresses  at 
different  points  on  the  body.  The  software  has  the  capability 
of  both  static  and  dynamic  analysis  of  a  structure,  although 
it  limits  the  number  of  degrees  of  freedom,  and  hence  the 
number  of  node  points,  for  certain  dynamic  analysis 
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techniques.  In  addition  to  tabular  output  of  nodal  locations, 
element  connectivity  and  nodal  deformations,  the  software  also 
has  the  capability  to  produce  graphical  output  for  analyzing 
model  geometry,  structural  deformations,  displacements  as 
functions  of  time  and  frequency,  animated  deformation  plots, 
and  stress  contour  plots. 

MSC/PAL2  requires  an  80386  (or  80486)  personal 
computer  with  at  least  a  20  megabyte  (Mb)  hard  disk  and  at 
least  2  Mb  of  random  access  memory,  with  additional  hard  disk 
space  being  required  for  very  large  problems .  Other 
requirements  include  at  least  one  floppy  disk  drive,  PC  DOS 
3.1  (or  higher)  operating  system,  an  80387  math  coprocessor 
and  either  a  CGA,  EGA,  or  VGA  monitor. 

2  .  Basic  Program  Operatior. 

MSC/PAL2  Version  4.0  is  a  menu  driven  program, 
beginning  at  the  top  level  with  six  main  menus.  Each  of  these 
menus  then  branches  out  into  varying  numbers  of  submenus, 
depending  on  the  application  chosen.  The  interrelationship  of 
these  menus  is  shown  in  Figure  4,  and  is  addressed  in  the 
discussion  that  follows. 

a.  File 

The  file  menu  has  only  two  choices.  The  first 
selection,  "Configure",  allows  the  user  to  change  the  disk 
drive  destination  to  which  the  database  files  will  be  saved. 
It  is  primarily  used  to  store  databases  on  alternate  drives 


12 


Figure  4  PAL2  Program  Execution 

(for  example  A,  B,  D)  for  very  large  problems.  The  second 
selection,  "Quit",  is  the  command  used  to  exit  the  program. 

b.  Edit 

The  Edit  menu  likewise  has  only  two  choices.  The 
first  choice,  "Text  file",  allows  the  user  to  view,  create  or 
edit  an  existing  text  file  utilizing  the  software's  text 
editor.  This  is  the  option  used  to  create  and  edit  the  model 
file  and  the  load  fil6(s)  which  are  utilized  in  the  analysis. 
The  second  choice,  "MSC-MOD",  is  an  interactive  graphical 
preprocessor  for  the  MSC,/PA1,2  program,  enabling  the  user  to 
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build  and  edit  two-dimeniiional  and  three-dimensional  finite 
element  models  representing  actual  structures.  Instead  of 
wciting  software  code  in  a  text  file,  MSC-MOD  uses  a  menuing 
system  overlaid  on  a  graphical  display  which  has  the  advantage 
of  letting  the  usei'  observe  on  the  display  the  results  of  a 
pt"  rticular  command.  MSC-MOD  can  be  used  in  place  of  MSC/PAL2 
to  create  a  PAL2  i:'odel,  or  it  can  be  used  to  edit  an  existing 
PAL2  model . 

c.  Analysis 

The  Analysis  menu  serves  the  purpose  of  providing 
commands  for  both  building  the  model  from  the  previously 
defined  model  text  file,  and  also  performing  the  actual  static 
and/or  dynamic  analyses  using  the  previously  defined  load 
file(s) .  The  "Build  Model"  submenu  generates  a  mathematical 
model  of  the  structure  defined  in  the  model  file,  which 
consists  of  a  set  of  second-order  ordinary  differential 
equations  of  motion,  representing  the  mass  and  stiffness 
mathematical  expressions  for  the  structure.  These  generated 
system  equations,  when  solved,  give  the  structure ' .s  nodal 
displacements  and  element  stresses.  Since  these  equations  are 
formed  from  the  mass,  stiffness  and  damping  of  each  element 
making  up  the  structure,  and  are  specified  for  each  nodal  DOF 
in  the  system,  matrices  are  employed  to  help  keep  track  of  all 
the  different  equations.  The  system  stiffness  matrix  is  built 
up  from  the  stiffnesses  of  each  element,  which  are  dependent 
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on  the  size,  shape,  tyj’e  and  material  make-up  of  each  element  . 
The  element  stiffnesses  are  defined  via  the  PAL2  "BEAM  TYPE" 
command,  fi.om  which  the  beam  cross-sectional  area,  bending 
moments  of  inertia,  aiid  toisional  moment  of  inertia  aie 
ccmi^uted.  The  system  mass  and  damping  matrices  are  computed 
by  the  program  in  an  analogous  manner,  with  the  mass  matrix 
originating  from  the  mass  density  specification  in  the  PAL2 
"MATERIAL  PROPERTIES"  conuuand,  and  the  damping  matrix 
originating  from  the  PAL2  commands  "DAMPER"  and  "DAMPING 
ELEMENT" . [Ref .  5: Section  1 , pp .  2-27  through  2-30] 

The  second  submenu  under  the  Analysis  heading  is 
the  "Static  Aiialysis"  menu  option.  Static  analysis  is  used  to 
determine  the  stiffness  and  element  stresses  of  the  structure 
when  subjected  to  some  static  load,  which  is  defined  in  a 
static  load  file.  The  static  loads  are  applied  to  the  model 
(built  using  the  "Build  Model"  submenu  command)  and  the 
program  computes  the  displacements.  Additional  loads  can  be 
applied  to  the  same  model,  if  desired,  without  building  the 
model  again.  In  fact,  the  analysis  can  be  run  repeatedly 
without  rebuilding  the  model.  After  all  of  the  loading  cases 
of  interest  have  been  run,  the  "Data  Recovery"  option  is 
invoked  in  order  to  obtain  element  stresses. [Ref.  6; Section  2, 
p.  3-8] 

The  third  submenu  under  the  Analysis  heading  is  the 
"Dynamic  .Analysis"  menu  option.  Dynamic  analysis  proceeds  in 
a  similar  manner  to  static  analysis,  although  multiple  load 
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conditions  cannot  be  run  on  the  same  model.  The  "Dynamic 
Analysis"  menu  option  is  very  flexible  in  that  it  allows  the 
user  to  compute  a  structure's  natural  frequencies  and  mode 
shapes  (using  the  "Normal  Modes  Analysis"  option) ,  responses 
to  frequency— dependent  forces  (using  the  "Frequency  Response 
Analysis"  option),  and/or  responses  to  time-varying  loads 
(using  the  "Transient  Response  Analysis"  option) .  Dynamic 
analysis  of  a  structure  should  be  run  only  after  the  static 
analysis,  and  the  normal  modes  analysis  should  bo  the  first 
dynamic  analysis  run  to  assess  t.he  overall  dynamic 
characteristics  of  the  structure .[ Ref .  5:Section  l,p.  3-4] 
d.  Graphics 

The  Graphics  menu  of  MSC/PAL2  has  two  submenus  from 
■which  to  choose:  "View"  and  "XYPlot"  .  The  "View"  option  can 
be  used  in  two  ways.  First,  it  can  be  used  to  check  model 
geometry,  prior  to  actually  analyzing  the  model,  to  make 
certain  that  all  node  points  are  properly  connected  with 
elements.  This  is  accomplished  using  the  "F2 :  Get  Dataset" 
option  from  the  first-level  "View"  menu.  "View"  is  also  used 
after  a  model  has  been  built  to  view  the  structural 
deformation  to  ^ome  load,  including  animation  of  the 
deformation.  This  is  accomplished  using  the  "FI:  Get  Model" 
option.  The  "View"  option  is  used  to  view  the  graphical 
results  of  both  static  and  dynamic  analyses,  in  both  the 
deformed  and  undeformed  cases. 
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The  second  option  under  the  Graphics  heading  is  the 
"XYPlot"  option.  This  option  is  primarily  used  to  graphically 
represent  frequency  and  transient  response  results  in  X-Y  plot 
format,  and  provides  graphical  displays  of  response  as  a 
function  of  time,  frequency,  element  number  and  subcase.  Like 
the  "View"  option,  "XYPlot"  is  menu-driven  and  uses  the  ten 
function  keys  to  select  the  different  plot  options. 

e.  Composite  Cap&bilities 

The  Composite  menu  is  a  pre-  and  postprocessor  to 
the  Analys.i.s  section,  used  to  cr  iate  composite  laminate 
descriptions  for  plate  material  properties  and  for  printing  or 
plotting  composite  results.  Laminates  with  up  to  200  layers 
made  from  up  to  ten  different  materials  can  be  made.  Since 
the  AIM-9  missile  is  not  of  composite  construction,  this 
section  of  the  MSC/PAL2  software  was  rot  used  in  this 
investigation.  It  is  mentioned  here  simply  for  completeness. 
[Ref.  6;Section  3,p.  1-1] 

f.  Advanced  CapahiJ-ities 

The  Advanced  menu  presents  some  additional 
capabilities  of  the  MSC/PAIj2  software,  and  includes  the 
following.  The  "Status"  menu  option  shows  whether  the  DOFs  of 
the  current  model  are  active,  eliminated  or  zeroed.  The 
"Equations"  menu  option  prints  out  the  terms  of  the  mass, 
stiffness  and  damping  matrices  for  the  most  recent  model 
created  by  the  "Build  Model"  option.  The  "Expanded  File"  menu 
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option  provides  expanded  MSC/PAIj2  model  data  sets,  which  are 
significantly  larger  than  the  original  data  sets  from  which 
they  were  created,  since  all  nodal  points  and  element 
connectivity  are  explicitly  specified.  The  "MSC/NASTRAN 
CONVERSION"  menu  option  translates  an  MSC/PA1j2  model  dataset 
to  MSC/NASTRAN  format.  Finally,  the  "REPLAY2"  menu  option, 
which  is  a  separate  program  from  MSC/P7iL2,  allows  the  user  to 
continuously  replay  text,  picture  and  animation  files  in  a 
slide-show  format .  Since  the  MSC/PAL2  program  was  only  used 
in  this  investigation  to  build  a  model  file  and  obtain  natural 
frequencies  and  mode  shapes,  these  advanced  capabilities  were 
not  utilized,  and  are  mentioned  here  simply  for  comp].eteness  . 
[Ref.  6; Section  4,pp.  1-1  through  6-1 j 
3 .  MSC/PAL2  Analysis  Procedure 

The  an?> lysis  of  a  generic  structure  using  the  MSC/PAL2 
finite  element  software  program  can  be  broken  down  into  three 
general  steps.  The  first  step,  system  definition,  is  the 
process  in  which  the  actual  structure  is  simulated  by  a 
mathematical  model,  defined  by  node  points  connected  with 
elements.  The  analyst  chooses  a  coord.- nate  system,  and  the 
locations  of  all  node  points  are  specified  relative  to  the 
origin  of  the  cocjrdinate  system  chosen .  After  all  of  the 
nodal  point  locations  are  defined  using  the  "NODAL  POINT 
LOCATIONS"  commands,  the  structure's  material  properties 
(Young's  modulus,  shear  modulus,  mass  density,  Poisson's  ratio 
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and  tensile  yield  stress)  are  specified  using  the  "MATERIAL 
PROPERTIES"  command.  Once  the  material  properties  have  been 
specified,  the  element  types  are  designated.  Element  type 
choices  include  several  types  of  curved  and  straight  beam 
elements,  quadrilateral  and  triangular  plate  elements,  and 
hexahedral,  pentahedral  or  tetrahedral  solid  elements. 
Following  the  element  type  specification  is  the  element 
connectivity  process,  using  the  "CONNECT",  "DO  CONNECT"  and/or 
"GENERATE  COtJNECTS"  commands  to  specify  the  nodal  points  that 
are  joined  by  each  element  in  the  model.  Finally,  the 
boundary  conditions  (for  example,  attachment  points)  are 
specified  using  the  "ZERO",  "ELIMINATE",  "ATTACH  AT"  and 
"OFFSET  CONNECTION"  commands .  After  the  model  has  been  fully 
defined,  the  "Build  Model"  menu  option  is  engaged  to  assenibla 
the  system  equations,  at  which  time  the  model  is  ready  for 
static  and/o.r  dynamic  analysis. 

The  second  general  step,  analysis,  is  the  process  in 
which  different  loading  conditions  are  applied  and  the 
structure's  static  and/or  dynamic  behavior  is  computed. 
Static  and  dynamic  loads  are  applied  to  the  already-built 
model  via  separate  load  files,  using  the  commands 
"DISELACEMENTS  APPLIED",  "ACCELERATIONS  APPLIED",  and  "FORCES 
AND  MOMENTS  APPLIED",  among  others,  and  the  resulting  nodal 
displacements  and  element  stresses  are  computed.  The  analysi.s 
is  conducted  using  the  "Static  Analysis"  and  "Dynamic 
Analysis"  submenus,  de.acribed  previously. 
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The  third  general  step,  graphical  postprocessing,  is 
the  process  in  which  the  static  and  dynamic  results  are 
presented  in  graphical  format.  The  "Graphics/View"  option  can 
be  used  to  show  the  structure's  undeformed  geometry, 
structural  deformation  and  deformed  animation.  The 
"Graphics/XYPlot"  option  can  be  used  to  plot  nodal 
displacements  and  accelerations,  and  also  element  stresses,  as 
a  function  of  time  or  frequency.  These  options  were  also 
previously  described  in  detail. 

For  the  conduct  of  an  actual  finite  element  analysis, 
these  three  general  steps  can  be  further  broken  down  into  a 
more  specific  plan  of  attack.  This  breakdown,  shown  in  Figure 
5,  consists  of  a  ten-step  process.  If  the  analyst  is 
conducting  a  full-blown  analysis,  complete  with  frequencies, 
mode  shapes,  stresses  and  deflections  caused  by  a  variety  of 
different  loading  conditions,  then  the  process  should  be 
completed  in  its  entirety,  if  possible.  If,  however,  the 
analyst  is  interested  in  only  certain  parameters  (for  example, 
natural  frequencies  and  mode  shapes) ,  then  only  those  steps 
relevant  to  that  particular  analysis  need  be  completed. 

4 .  MSC/PAli2  liiiaitations 

The  maximum  problem  size  for  this  version  of  PAL2  is 
dependent  upon  how  many  DOFs  are  in  the  model,  as  well  as  the 
type  of  DOF.  Initially,  there  are  G  global  DOFs  in  each 
model,  where  G  is  six  times  the  number  of  nodal  points,  since 
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Figure  5  MSC/PAL2  Suggested  Order  of  Analysis 


each  nodal  point's  motion  can  be  represented  by  three 
translations  and  three  rotations  about  the  coordinate  axis 
system  chosen.  The  global  DOFs  are  partitioned  into  the 
zeroed  DOFs  (created  with  the  "Build  Model"  command  "ZERO") 
and  the  free  DOFs.  The  free  DOFs  are  further  partitioned  into 
the  eliminated  DOFs  (created  with  the  "Build  Model"  command 
"Eliminate")  and  the  active  DOFs.  With  this  in  mind,  the 
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problem  size  limitations  for  this  version  of  the  MSC/PAL2 
software  are  as  follows  [Ref.  5:p.  A-2] : 


•  Build  Model:  2000  nodes,  12000  DOF 

•  Statics:  12000  free  DOF  or  225  active  DOF 

•  Dynamics :  Normal  modes 

Subspace  iteration:  12000  active  DOF 
QR  method:  225  active  DOF 
Jacobi,  method:  225  active  DOF 
Transient  response;  225  active  DOF 
Frequency  response;  150  active  DOF 

If  a  model  has  more  active  DOFs  than  shown  here,  and  the 
analyst  wishes  to  conduct  a  detailed  analysis,  then  the  number 
of  active  DOFs  must  be  reduced  r'ia  static  condensation  (for 
static  analysis)  or  Guyan  reduction  (for  dynamic  analysis) . 
MacNeal-Schwendler  [Ref.  5; Section  l,p.  2-43]  addresses  in 
some  detail  the  procedure  for  accomplishing  this  reduction. 
Since  the  main  thrust  of  developing  the  AIM-9  model  was  to 
analytically  determine  the  ir.issile's  natural  frequencies  and 
mode  shapes,  the  normal  modes  subspace  iteration  method  was 
used,  and  no  reduction  was  necessary. 

C.  MSC-PAL2  SIDSWINDSR  MISSILE  AHALY7ZCAL  MODEL 
1 .  Model  Geometry 

The  MSC/PAL2  software  code  written  to  define  the  AIM-9 
geometry  is  shown  in  Appendix  A.  The  actual  missile  geometry 


resulting  from  this  code  is  shown  in  Figure  6,  both  with  the 
hidden  elements  showing  and  with  them  not  showing.  From  the 
code,  it  can  be  seen  that  a  total  of  706  node  points  w»are  used 
to  define  this  structure;  610  node  points  for  the  missile 
body,  36  for  the  front  fins,  and  60  for  the  rear  wings.  This 
large  number  of  points  was  required  to  describe  the  missile 


Figure  G  MSC/PAL2  AIM-9  Geometry  Model 


for  several  reasons.  fcirst,  instead  of  using  curved  beam 
segments  (which  cannot  be  graphically  displayed  in  the  PAL2 
.software)  to  make  up  the  cylindrical  missile  tody,  concentric 
rings  of  16  nodes  joined  with  small  quadrilateral  plate 
elements  were  spaced  three  inches  apart  to  define  the  missile 
body.  Secondly,  these  rings  were  spaced  this  closely  because 
of  the  recommendation  by  MSC  [Ref.  7]  to  avoid  aspect  ratios 
(length/width  ratios)  in  the  plate  elements  larger  than  four 
to  one.  Thirdly,  a  fairly  uniform  spatial  distribution  of 
node  points  is  required  to  properly  compute  resonant 
frequencies  and  mode  shapes.  Since  the  missile  model  was 
built  using  the  actual  dimensions  of  the  missile,  tne  result 
was  35  concentric  rings  of  node  points  making  up  the  uniform 
cross-section  of  the  missile  body  and  three  concentric  rings 
of  node  points  making  up  the  tapered  nose  cone,  all  tapering 
into  one  node  point  representing  the  nose  of  the  missile. 
Quadrilateral  plate  elements  of  0.25  inch  thickness  were  used 
to  define  the  missile  body  and  the  tapered  nosecone,  while 
triangular  plate  elements  of  0.25  inch  thickness  were  used  for 
the  tip  of  the  nose  cone  and  the  back  plate.  The  three  points 
where  the  ATM-9  missile  attaches  to  the  LAU-7  launcher  rail 
were  modeled  by  zeroing  the  displacements  at  the  applicable 
node  points.  Translation  in  the  x,  y,  and  z  directions  and 
rotations  about  the  y  and  z  axes  were  zeroed  for  node  points 
85,  213,  and  373.  Rotation  about  the  x  axis  was  not  zeroed. 
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since  motion  in  the  pitch  direction  is  the  primary  motion  ,of, 
interest . 

After  the  missile  body  was  constructed,  the  front  fins 
were  added  using  a  combination  of  quadrilateral  and  triangular 
plates,  each  of  0.10  inch  thickness,  with  great  care  taken  to 
keep  the  individual  plate  aspect  ratios  at  less  than  four  to 
one.  Connecting  the  front  fins  to  the  missile  body  proved  to 
be  a  problem  because  of  the  fact  that  the  AIM-9' s  front  fins 
are  actually  movable  canards  which  exhibit  rotation  about 
their  attactiment  points  to  provide  missile  guidance.  The  PAL2 
software  has  no  provisions  for  modeling  this  type  of 
arrangement,  so  the  fins  were  attached  to  the  missile  body  via 
a  small  0.35  inch  plate,  thereby  prohibiting  the  rotation  of 
the  fins  in  the  x-y  plane  about  their  attachment  points.  This 
resulted  in  a  stiffness  problem,  with  the  fins  rotating  about 
the  2-axis.  Additional  work  is  required  in  this  area  to 
reflect  a  more  representative  attachment  method. 

The  rear  wings  were  the  next  building  block  added  to 
the  model.  These  wings  were  modeled  using  a  combination  of 
0.5  inch  thick  quadrilateral  and  triangular  plates,  and  were 
rigidly  attached  to  the  missile  body  using  previously  defined 
missile  body  node  points. 

All  that  remained  to  complete  the  missile  model  was 
the  addition  of  mass  parameters  to  represent  the  motor, 
warhead,  target  detector  and  guidance  and  control  components 
of  the  AIM-9  missile.  This  was  accomplished  using  point 


masses  located  at  the  midpoint  of  each  component's  length, 
with  the  exception  of  the  motor  component,  which  was 
represented  by  two  point  masses  because  of  its  70  inch  length. 

Finally,  the  material  coinposit ron  of  each  element  in 
the  model  was  defined  using  generic  steel  values  of  Young's 
modulus,  shear  modulus,  mass  density,  Poisson's  ratio  and 
tensile  yield  strength  for  the  missile  body  and  front  fins, 
and  generic  aluminum  values  for  the  rear  wings .  These 
material  properites,  along  with  the  nodal  locations  and 
element  connectivity,  are  used  to  generate  equations  that, 
when  solved,  give  the  system  displacements  and  stresses. 

The  completed  AIM-9  missile  model  consisted  of  706 
nodal  point  locations,  which  resulted  in  a  total  of  4221  DOFs. 
The  model  was  constructed  of  717  elements  of  three  different 
types;  660  quadrilateral  plate  elements,  52  triangular  plate 
elements  and  five  point  mass  elements.  Quadrilateral  plate 
elements  were  used  as  much  as  possible  because  they  are  less 
stiff  than  the  triangular  plate  elements,  thereby  giving 
better  and  more  accurate  results  [Ref,  5: Section  2,p.  2-3]. 

2 .  Determination  of  Modal  Parameters 

Upon  completion  of  the  model  geometry,  and 
verification  of  the  nodal  point  locations  and  connectivity 
using  the  "Graphics-View"  option,  the  "Build  Model"  menu 
option  was  invoked  to  form  the  system  stiffness,  damping  and 
mass  matrices.  These  matrices  arise  from  Newton's  Second  Law 


of  Motion,  upon  which  all  vibration  work  is  based.  Written  as 
a  second-order  differential  equation  in  the  time  domain, 
Newton's  Second  Law  can  be  written  as  follows; 

Mx(t)  +  Cx<t)  +  Kx(t)  =  f(t)  (1) 

where  M  =  system  mass  maurix 

C  =  system  damping  matrix 
K  =  system  stiffness  miat  rix 
x(t)  =  nodal  displacement  vector 
x(t)  =  nodal  velocity  vector 
’x(t)  =  nodal  acceleration  vector 
f(t)  =  applied  force  vector 

After  the  modei  was  built,  a  riormal  modes  analysis  was 
conducted  to  determine  the  missile' s  natural  frequencies  and 
mode  shapes.  The  system  equations  formed  for  a  normal  modes 
analysis  include  the  system  stiffness  and  mass  matrices, 
written  in  matrix  notation  as  follows: 

Mx  +  Kx  =  0  (2) 

The  right-hand  side  of  this  equation  is  set  equal  to  zero 
since  there  is  no  forcing  in  normal  modes  analysis.  All 
motion  is  assumed  to  be  steady— state  harmonic  motion.  The 
structure,  in  this  case  the  AIM— 9  missiJ.e,  vibrates  at  some 
frequency,  f,  with  each  point  moving  either  in  phase  with  the 
other  points  or  out  of  phase.  Therefore,  the  nodal 
acceleration  vector  may  be  written  in  terms  of  the 
displacement  vector  as  follows: 

X  =  -(27rf)^x  (3) 
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Substituting  this  into  Equation  2  gives: 

[K  -  (27Uf)='M]x  =  0  (4) 
In  the  case  of  the  AIM-9  missile,  the  dimensions  of  the  mass 
and  stiffness  matrices  is  4221  x  4221,  since  the  model  has 
4221  active  DOFs.  One  can  immediately  see  the  necessity  of 
using  a  computer  program  for  this  type  of  analysis,  as  the 
manual  manipulation  of  this  many  equations  would  be  virtually 
impossible.  The  nontrivial  solution  to  Equation  4  yields  the 
system  resonant  frequencies,  each  of  which  has  associated  with 
it  a  chai-acteristic  vibration,  or  mode  shape.  Since  there  is 
no  forcing  in  a  normal  modes  analysis,  the  scaling  values  for 
each  mode  shape  are  arbitrary,  and  only  the  relative  shape  is 
important . [Ref .  5:Section  l,p.  2-40] 

A  normal  modes  analysis  can  be  accomplished  by  three 
different  methods.  The  default  method  of  eigenvalue 
extraction,  the  QR  method,  is  much  faster  than  the  Jacobi 
method,  which  is  a  more  accurate  method  for  unrestrained 
structures .  Both  of  these  methods  require  the  system 
equations  to  be  reduced  in  size  to  the  point  where  there  are 
no  more  than  225  active  DOFs.  This  reduction  is  accomplished 
with  the  PAL2  commands  "Eliminate"  and  "Activate" .  Since 
typically  only  the  lowest  few  modes  are  excited  to  any 
appreciable  degree  by  dynamic  forces,  these  are  the  modes  of 
primary  interest  to  the  designer  or  analyst.  The  subspace 
iteration  method  is  the  method  of  choice  when  only  the  lowest 
frequencies  are  needed,  because  it  accomodates  up  to  12000 
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active  DOFs  and  no  reduction  ot  the  system  equations  is 
required.  This  method  solves  for  the  eigenvalues  and 
eigenvectors  with  the  mass  and  stiffness  matrices  in  the  form 
of  Equation  4,  obtaining  these  equations  from  the  physical 
properites  of  the  structure  coded  into  the  model,  and  was  the 
method  used  for  conducting  a  normal  modes  analysis  of  the  AIM- 
9  missile  model.  The  results  of  this  analysis  are  given  in 
Appendix  B . 

3 .  Modification  of  Modal  Parameters 

The  modal  parameters  of  the  AIM-9  model,  specifically 
the  natural  frequencies  and  mode  shapes,  may  be  readily 
modified  to  reflect  alternative  arrangements  of  the  missile 
configuration,  simply  by  altering  the  mass,  stiffness  ox‘ 
damping  properties  input  into  the  model.  In  the  construction 
of  the  missile  model,  various  ditierent  configurations  were 
analyzed  to  determine  the  effect  of  changing  these  properties. 
The  results  of  the  five  cases  analyzed  are  presented  and 
discussed  in  Appendix  C.  It  can  be  seen  from  these  results 
that  the  addition  or  movement  of  the  point  masses  can  have  a 
substantial  effect  on  both  the  natural  frecjuencies  and  the 
mode  shapes.  Adding  mass  to  the  structure  shifts  the 
resonance  to  a  lower  frequency,  as  does  decreasing  the 
stiffness,  while  removing  mass  shifts  the  resonance  to  a 
higher  frequency,  as  does  increasing  the  stiffness.  This  is 
helpful  in  avoiding  a  known  excitation.  Moving  the  location 
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ot  a  point  mass  only  a  shoit  distance  had  the  effect  of 
completely  changing  the  mode  chapes .  Chaiiging  the  material 
properties  of  parts  of  the  structure  also  affected  the 
response.  Since  these  results  were  obtained  from  a  normal 
modes  analysis  in  which  only  the  relative  shape  of  the  mode  is 
important  and  not  the  amplitude,  damping  is  not  addressed. 
However,  damping  would  need  to  be  increased  to  reduce  the 
response  at.  a  j'art icul ar  re.sonance  . 

Knowing  the  effects  of  changing  the  various  properties 
of  a  model  enables  the  designer  or  analyst  to  investigate 
different  design  iterations  or  feasibility  studies  prior  to 
making  any  hardware  changes.  For  example,  if  ascertaining 
the  effect  of  adding  more  fuel  to  the  AIM-9  missile  to  extend 
its  range  was  desired,  mass  would  be  added  to  the  existing 
model  to  represent  the  increa.sed  fuel  required,  and  additional 
node  points  would  be  added  to  represent  the  increased  length 
of  the  fuel  tank  component.  A  new  normal  modes  analysis  could 
then  be  run,  and  the  new  natural  frequencies  and  mode  shapes 
obtained.  Another-  possible  missile  upgrade  includes 
incorporating  a  new  type  of  seeker  head.  The  effects  of  this 
potential  modification  could  also  be  observed  quite  easily. 
For  these  studies  to  be  of  value,  it  is  imperative  that  the 
baseline  model  accurately  represent  the  actual  structure  not 
just  geometrically,  but  in  the  modal  responses  as  well.  For 
this  reason,  it  is  important  to  utilise  the  combined  testing 
and  analysis  approach  to  fix'st  obtain  the  missile's  modal 
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parameters  experimentally,  then  modify  the  analytical  model, to. 
match  the  experimental  results.  It  is  at  this  point  that  the 
design  iterations  or  feasibility  studies  can  be  conducted  on 
the  analytical  model  with  a  reasonable  amount  of  confidence  in 
the  results  obtained. 
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III.  EXPERIMENTAL  MODAL  ANALYSIS 


A.  GENERAL 

III  the  same  vein  that  the  finite  element  method  can  be 
considered  a  "forward-type"  analysis,  in  that  a  structure's 
constituent  properties  are  input  and  the  response  is 
predicted,  modal  analysis  can  be  considered  an  "inverse-type" 
analysis,  in  that  a  structure's  response  to  an  excitation  is 
measured,  and  from  the  re.sulting  transfer  function,  the 
structure' s  modal  parameters  can  be  obtained  using  analytical 
parameter  estimation  methods.  The  .".ajor  advantage  of 
conducting  a  modal  analysis  of  a  structure  is  that  simplified 
mass-spring-damper  models  can  be  constructed  with  a  good 
degree  of  certainty  from  the  actual  measured  response  of  the 
structure.  The  drawback  of  modal  testing  is  that  design 
iterations  or  feasibility  studies  cannot  be  accomplished 
unless  the  structure  is  physically  modified.  The  best 
approach  is  to  combine  the  two  methodologies,  by  conducting  a 
modal  analysis  of  the  test  structure  to  obtain  the  actual 
response,  and  then  use  this  information  to  make  the  PAL2  model 
more  accurate  for  use  in  future  design  iterations.  This 
combined  testing  and  analysis  approach  is  shown  in  block 
diagram  form  in  Figure  7  [Ref.  2:p.  3-4], 


MODAL  FINITE 


Figure  7  Combined  Testing  and  Analysis 

Structural  Measurement  Systems,  Inc.  [Ref.  2:p.  1-2] 
defines  modal  analysis  as  the  process  of  characterizing  the 
dynamic  properties  of  an  elastic  structure  in  terms  of  its 
modes  of  vibration.  In  the  majority  of  present  day  modal 
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analysis  practice,  as  well  as  in  finite  element  analysis  as 
mentioned  previously,  the  motion  of  a  physical  system  can  be 
assumed  to  be  adequately  described  by  a  set  of  simultaneous 
second-order  linear  differential  equations  of  the  form: 

Mx(t)  +  Cx(t)  +  Kx(t)  =  f(t)  (5) 

If  a  system  has  n-DOFs,  then  the  above  vectors  are  n- 
dimensional  and  the  matrices  are  nxn .  These  system  equations 
are  transformed  into  their  equivalent  transfer  function  form 
by  taking  the  Laplace  transform,  resulting  in  the  relation: 

X(s)  =  H(s)  F(s)  (6) 

where  H(s)  is  cal''ed  the  transfer  matrix,  which  can  then  be 
written  in  terms  of  modal  frequency,  damping  and  modal 
vectors.  It  is  important  to  note  that  each  row  and  column  of 
the  transfer  function  matrix  contains  the  same  modal  vector, 
multiplied  by  some  component  of  itself.  Therefore,  only  one 
row  or  column  of  the  transfer  matrix  needs  to  be  measured  in 
order  to  identify  all  the  modal  parameters  of  a  structure, 
provided  the  following  assumptions  are  met .  The  first 
assumption  is  that  the  motion  of  the  structure  can  be 
described  by  linear  second-order  equations.  The  second 
assumption  is  that  symmetry  of  motion  exists,  described  by 
Maxwell's  Reciprocity  Theorem,  which  implies  that  the  FRF 
measured  between  any  two  DOFs  is  independent  of  which  of  them 
is  used  for  excitation  or  response.  The  third  and  final 
assumption  is  that  no  more  than  one  mode  exists  at  each  pole 
location  of  the  system  transfer  matrix. [Ref.  8] 
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Modal  analysis  is  conducted  on  a  structure  for  various 
reasons .  Perhaps  one  of  the  most  important  is  to  verify  and 
help  improve  the  accuracy  of  finite  element  models. 
Experimental  measurements  obtained  on  the  structure  being 
analyzed  result  in  a  physical  check  of  the  accuracy  of  the 
mathematical  model.  If  the  model  predicts  the  same  modes  and 
frequencies  as  those  obtained  experimentally,  then  the  model 
is  considered  to  be  accurate.  Otherwise,  the  model  must  be 
modified  to  match  the  same  response  as  that  obtained 
experimentally,  assuming  the  experimental  response  is  correct, 
other  important  reasons  for  doing  modal  testing  and  analysis 
include  troubleshooting  noise  and  vibration  problems,  quickly 
evaluating  fixes  made  to  structures  in  order  to  solve  noise 
and  vibration  problems,  and  formulating  dynamic  models  for 
parts  of  structures  that  are  too  difficult  or  time-consuming 
to  model  analytically.  In  this  investigation,  a  modal 
analysis  test  setup  was  completed  to  be  used  to  analyze  the 
response  of  the  structure  of  interest  (AIM-9  missile)  in 
conjunction  with  another  structural  member  (P—3  outer  wing) . 

B.  PREVIOUS  WORK  ACCOMPLISHED 
1 .  Test  Setup 

Work  was  begun  on  this  effort  by  Hollyer  [Ref.  1]  who 
designed  and  built  the  GVT  stand,  which  was  then  used  to 
conduct  a  vibration  characterisation  of  the  AIM-9  missile  (in 
pitch  only) ,  employing  forced  oscillatory  inputs  over  a 
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frequency  range  of  0-50  hz .  In  this  test,  Hollyer  used  analog 
instrumentation  (vice  a  digital  spectrum  analyzer)  in  which  a 
shaker  assembly  was  utilized  in  conjunction  with  a  function 
generator  to  excite  the  structure  sinusoidally.  To  measure 
the  missile's  response,  two  piezoelectric  accelerometers  were 
used,  one  mounted  rigidly  in  series  with  the  shaker  assembly 
(the  reference  transducer)  and  one  roving  over  the  missile  to 
read  magnitude  and  phase  shift  at  various  locations  on  the 
missile  (the  response  transducer)  using  an  oscilloscope. 

2 .  Previous  Results 

Using  the  above  apparatus,  the  missile  system  was 
excited  in  pitch,  with  the  excitation  force  varying  from  10  to 
40  pounds  peak-to-peak .  Pitch  modes  were  observed  over  a  band 
of  24.4  to  25.6  hz  and  near  34  hz,  and  a  heave  mode  was 
observed  near  38  hz .  The  first  elastic  mode  of  the  missile, 
obtained  by  placing  the  shaker  at  the  node  of  the  first  rigid 
mode,  appeared  to  occur  at  about  52.5  hz .  However,  this  mode 
was  terribly  polluted  with  the  motion  of  the  GVT  stand  itself. 


Exact  frequencies  could 

not  be  determined 

due  to 

the 

inaccuracies 

of  the  data 

collection 

equipment 

used . 

The 

results  from 

this  test  are  shown  in 

Figure  8  . 

Since 

the 

primary  thrust  of  this  effort  was  to  determine  the  natural 
frequencies  and  mode  shapes  of  the  missile  in  pitch,  the  GVT 
was  built  and  stiffened  for  excitation  in  the  longitudinal 
direction  only.  Even  so,  the  system  was  excited  in  sway  and 
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yaw,  with  a  sway  mode  observed  at  about  7.2  hz  and  a  yaw  mo'Je' 
at  about  9.8  hz  .  Hov;ever,  these  modes  are  very  suspect  due  to 
significant  motion  of  the  GVT  during  excitation.  To  obtain 
more  accurate  missile  response  in  sway  and  yaw,  the  GVT  would 
need  to  be  stiffened  laterally . [Ref .  1] 


POSITION  AFT  OF  THE  NOSE  OF  THE  MISILE 
/iwrwci;» 


Figure  8  Missile  Modes  for  Single-Rail  AIM-9  [Ref.  1] 


After  the  above  modes  were  obtained,  the  missile 
system  vras  modeled  as  the  simplified  spring  lumped  mass  model 
showii  in  Figure  9.  This  model  was  represented  as  a  two  lumped 
mass  system  possessing  no  internal  damping,  with  two  springs 
to  provide  the  two  degrees  of  freedom  to  model  the  lowest 
resonance  of  the  pitch  and  heave  modes.  The  theory,  logic  and 
thought  process  beliind  this  model  can  be  found  in  more  detail 
in  Hollyer' s  investigation  [Ref.  1] . 
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Figure  9  Spring  -  Mass  Model  [Ref.  1] 

C.  MODAL  TEST  SETUP  CONSIDERATIONS 
1 .  Transducer  Selection 

Prior  to  conducting  a  modal  test,  a  multitude  of 
cc^iSiderations  must  be  taken  into  account  -  Perhaps  one  of  the 
most  important  considerations  is  the  selection  of  a 
transducer.  A  transducer  is  a  device  which  converts  the 
vibratory  motion  of  the  structure  into  an  analogous  electrical 
signal.  Transducers  come  in  three  basic  types:  displacement 
*■ '  .  jd'!--  ,  velocity  transducers  and  accelerometers.  The 

accelerometer  is  by  far  the  most  common  and  is  replacing  the 
velocity  transducer  in  the  majority  of  today' s  applications 
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because  of  its  superior  phase  and  frequency  response.  One  of 
the  most  common  accelerometers  in  use  today  is  the 
piezoelectric  accelerometer,  which  is  a  small,  lightweight 
transducer,  useable  over  a  wide  frequency  range.  Since  the 
output  of  a  piezoelectric  accelerometer  is  a  very  low  voltage, 
a  charge  amplifier  must  be  used  to  amplify  the  signal  to  a 
high  enough  level  to  be  displayed  by  the  spectrum  analyzer. 
A  recent  development  in  accelerometers  is  the  integrated- 
circuit-piezoelectric  (ICP)  concept,  in  which  amplifiers  in 
the  form  of  microelectronic  integrated  circuits  are  built  into 
the  transducer.  This  precludes  the  need  for  a  charge 
amplifier,  but  a  constant  current  power  supply  is  required  to 
operate  an  ICP  transducer. 

Considerations  in  selecting  an  accelerometer  for  a 
specific  project  include  sensitivity,  frequency  range,  size, 


and  long-term 

stability  [Ref.  9] . 

If  one  or 

more  of  these 

properties 

is 

important  for  a 

particular 

application. 

tradeoffs 

may 

have  to  be  made. 

and  some 

of  the  other 

properties  may  have  to  be  sacrificed.  Accelerometers  are 
usually  selected  to  have  a  sufficiently  high  resonance 
frequency  so  that  the  frequency  response  can  be  considered 
flat  throughout  the  operating  frequency  range.  Accordingly, 
the  accelerometer  reference  frequency  should  be  at  least  five 
times  the  maximum  frequency  of  interest  [Ref.  10]  .  By 
following  this  rule,  slight  variations  in  resonance  frequency 
due  to  different  mounting  methods  can  be  ignored. 


2.  Transducer  Mounting 

The  primary  consideration  in  selecting  a  transducer 
mounting  method  is  that  tha  mounting  method  does  not  introduce 
any  impedance  into  the  system  that  would  result  in  errors  in 
the  generated  data.  Various  mounting  methods  are  available, 
including  stud  mounting,  wax  mounting,  cementing  and  hand-held 
probes  and  magnets.  The  preferred  method  is  to  use  an 
isolated  stud  mount  whereby  the  transducer  is  physically 
attached  to  the  structure  via  a  screw  stud.  This  method 
provides  the  best  frequency  range,  accuracy  and  long-term 
stability,  and  problems  with  dust,  moisture  or  temperature  are 
normally  not  encountered.  However,  this  is  a  somewhat 
permanent  method  which  modifies  the  test  structure.  If  the 
test  structure  is  one  in  which  drilled  and  tapped  holes  are 
not  desired  (for  example,  a  missile)  ,  the  next  best 
alternative  is  beeswax.  This  miethod  allows  for  quick  setup 
and  transducer  movement,  while  still  providing  for  a  moderate 
frequency  range  and  accuracy.  Manufacturer  testing  [Ref.  11] 
has  shown  that  frequency  response  equivalent  to  stud  mounting 
can  be  achieved  with  a  thin  coating  of  wax  on  clean,  flat 
surfaces.  If  the  surface  is  curved,  a  composite  mounting  base 
can  be  used,  which  can  easily  be  sanded  for  mounting  to 
contoured  surfaces.  Cementing  is  sometimes  difficult  because 
of  the  setting  time  required,  and  also  because  of  adhesion 
problems  to  certain  surfaces.  Magnets  allow  for  quick  setup 
and  movement,  but  this  method  is  frequency-limited  due  to  the 
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lack  of  a  rigid  attachment.  Al.so/  mass  loading  of  small 
structures  may  occur.  Hand-held  probes,  while  being  the 
simplest  method,  are  not  recommended. 

3 .  Transducer  Calibration 

It  is  not  necessary  to  calibrate  the  transducers  used 
in  a  modal  test  if  all  that  is  desired  is  the  acquiring  of  a 
structure's  mode  shapes.  However,  transducer  calibration  is 
essential  when  developing  a  modal  model  to  be  used  for 
structural  modification  studies.  There  are  basically  three 
types  of  transducer  calibration  that  can  be  used.  Absolute 
calibration  of  a  single  transducer  requires  a  known  source 
(for  example,  Ig  at  IKhz)  which  is  fed  through  the  transducer 
and  read  directly  on  a  spectrum  analyzer.  Referential 
calibration  requires  a  well-calibrated  transducer  to  be  used 
as  a  standard  to  which  the  user's  transducers  are  referenced. 
Finally,  the  constitutive  or  ratio  calibration,  which  combines 
aspects  of  the  absolute  and  referential  methods,  uses  some 
constitutive  law  (for  example,  dividing  all  measurements  by 
the  ratio  (A^oit./fvoit.)  *  m)  to  relate  a  pair  of  transducers  to 
each  other. [Ref.  2:p,  10-2] 

4 .  Excitation  Techniques 
a.  Shaksr  Excitation 

The  objective  of  an  exciter  in  a  modal  test  is  to 
input  a  force  into  the  structure  under  test  to  excite  all  the 
modes  of  interest .  VariovJS  excitation  techniques  have  been 
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developed  to  help  the  analyst  solve  a  particular  vibration 
problem.  One  of  the  most  common  techniques  is  to  use  a 
hydraulic  or  electro-mechanical  shaker  assembly.  A  typical 
shaker  setup  for  a  modal  test  is  shown  in  Figure  10.  The 
shaker  should  be  attached  to  the  test  structure  in  such  a  way 
that  excitation  of  the  structure  occurs  at  a  single  point  and 
in  a  single  direction.  A  thin  flexible  "stinger"  with  a  force 
transducer  mounted  at  the  end  can  be  used  to  reduce  bending 
moments .  Using  a  shaker  assembly  requires  a  certain  amount  of 
fixturing,  but  it  gives  good  control  of  the  excitation 
bandwidth  and  level,  and  allows  several  choices  of  excitation 
signals.  These  signals  are  of  two  primary  types: 
deterministic  excitation  (which  includes  sinusoidal  sweep, 


^  ^ 
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excitation  (which  includes  pure  random,  pseudo  random, 
periodic  random  and  burst  random) .  Each  of  these  types  of 
excitation  techniques  ha.s  its  advantages  and  disadvantages. 
These  characteristics  will  not  be  listed  here,  but  are 
discussed  in  great  detail  by  SMS,  Inc.  [Ref,  2:pp.  9—2  through 
9-44]  .  The  point  is  that  the  analyst  must  select  an 
excitation  technique  based  on  the  considerations  and 
constraints  of  a  particular  application.  Examples  of  these 
considerations  are  test  measurement  time,  signal  leakage 
control,  bandwidth  control,  level  of  signal-to-noise  ratio  and 
distortion  removal . 
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Signal  Gcnctalor 


Figure  10  Typical  Shaker  Setup 

b.  Impulse  (Impact)  Excitation 

The  latest  advancement  in  medal  analysis  is  the  use 
of  an  impulse,  or  tap,  hanmier  to  impart  broadband  excitation 
to  the  test  structure,  thus  exciting  all  modes  simultaneously. 
The  excitation  force  response  can  be  altered  by  changing  the 
hammer  tip  material,  which  is  determined  by  the  frequency 
range  of  interest.  A  soft  rubber  tip  will  result  in  a  longer 


imiiact  time  which  will  concentrate  the  input  energy  into  the 
lower  frequency  ranges.  Conversely,  a  hard  steel  tip  will 
result  in  shorter  impact  times  for  cases  where  a  higher 
frequency  range  is  of  interest.  The  impact  technique  has 
several  advantages  over  a  shaker  assembly  in  that  it  is  fast, 
easy,  low  cost,  convenient  and  portable,  and  it  does  not 
require  any  external  signal  source.  However,  it  provides  only 
a  low  input  energy  with  no  randomization  of  the  input,  a  low 
signal-to-noise  ratio,  and  no  control  of  the  frequency  content 
if  a  particular  excitation  bandwidth  is  desired.  Corelli  and 
Brown  [Ref.  12]  state  that,  in  general,  impact  excitation  is 
the  worst  possible  excitation  method  to  use  for  most 
structural  tests,  as  the  technical  disadvantages  far  outweigh 
the  technical  advantages.  However,  the  convenience  and 
quickness  of  impact  testing  is  so  dominant  that  it  generally 
outweighs  the  technical  disadvantages  and  is  often  chosen  as 
the  excitation  technique. 

5.  Specifying  the  Number  of  Degrees -of -Freedom 

As  stated  earlier,  a  DOF  is  a  measurement  point,  and 
direction  defined  on  a  structure.  Any  free  point  on  a 
structure  generally  has  six  DOFs;  three  translational  and 
three  rotational.  When  deciding  oii  the  nun\ber  of  DOFs  to  use 
for  a  particular  test,  the  analyst  must  consider  the  purpose 
of  the  test,  and  the  geometry  of  the  actual  structure  under 
test.  If  the  purpose  of  the  modal  test  is  simply  to  verify 
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modal  frequencies  previousl;^  preaicted  analytically,  then  only 
a  few  DOFs  would  be  required.  However,  if  the  purpose  of  the 
test  is  to  construct  a  mathematical  model  interpolated  from 
modal  measurements  at  discrete  locations,  then  sufficient  DOFs 
must  be  used,  taking  the  structure's  geometry  into  account,  to 
ensure  that  all  mode  shapes  in  the  frequency  range  of  interest 
are  distinguishable  and  linearly  independent.  It  is  the 
geometrical  complexity  of  a  structure's  mode  shapes  which 
ultimately  determines  the  nuniber  of  DOFs  required.  Known  mode 
shapes  of  similar  structures  tested  previously  should  be  used 
whenever  pos.sible  to  determine  how  many  DOFs  with  which  to 
begin.  Care  should  be  taken  in  not  using  more  DOFs  than  are 
necessary  to  define  the  mode  shapes  for  a  structure,  as  this 
results  in  longer  computation  and  test  time. 

D  SIDEWINDER  MISSILE  MODAL  ANALYSIS  TEST  SE.'^UP 
1 .  Test  Stxructure 

The  GVT  stand  utilized  in  the  missile  modal  analysis 
setup  is  located  on  the  basement  floor  of  Kalligan  Hall  at  the 
Naval  Postgraduate  School  in  Nonterey,  California.  The  GVT  is 
the  same  structure  const.ucted  and  used  by  Hollyer  [Ref.  1], 
except  for  the  addition  of  a  dozen  additional  fasteners 
joining  the  top  and  bottom  plates  to  the  8  inch  by  8  inch  I 
beams,  and  two  14  inch  bolts  joining  the  top  and  bottom  plates 
tovjether  for  the  purpose  of  providing  further  stiffening  of 
the  structure.  This  structure  was  initially  designed  to 
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provide  minimal  notion  in  the  frequency  range  of  interest  (0 
to  50  hz)  to  longitudinal  test  member  excitation,  with  an 
initial  resonant  frequency  of  approximately  52,5  hs .  The 
structure  has  no  added  stiffening  for  excitation  in  the 
lateral  direction. 

2 .  Test  Equipment 

A  basic  modal  test  system  is  usually  comprised  of  a 
spectrum  analyzer,  a  computer  with  associated  modal  analysis 
software,  a  shaker  or  impact  hammer  to  provide  the  excitation 
(with  an  external  function  generator,  if  needed) ,  a  force 
transducer  to  provide  driving  point  measurements,  response 
transducers  to  measure  the  response  to  the  excitation,  and  the 
actual  structure  to  be  tested.  The  Sidewinder  missile  modal 
analysis  test  setup  is  depicted  in  block  diagram  form  in 
Figure  11,  and  the  actual  setup  is  shown  in  Figure  12. 
Specific  details  of  the  individual  components  of  this  test 
setup  are  discussed  in  Appendix  D. 

The  spectrum  analyzer  used  in  this  test  setup  is  the 
Scientif ic--Atlanta  Model  SD-380,  with  the  four  channel  opcion 
installed.  This  option  allows  the  parallel  sampling  of  four 
channels  of  simultaneous  inforraation,  thus  permitting  the 
force  transducer  to  be  connected  to  channel  A,  and  the  three 
Flexcel  response  transducers  to  be  connected  to  channels  B,  C, 
and  D.  The  SD-380  used  in  the  missile  test  setup  also  has  a 
Synchronous  Signal  Generator  (S3G)  capability,  precluding  tha 


Figuro  11  AIM-9  Missile  Modal  Analysis  Block.  Diagram 
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need  for  a  separate  function,  generator  to  send  the  excitation 
signal  to  the  shaker.  The  SSG  signal  output  types  available 
include  broadband  white  noise,  pseudo  random  noise,  band 
limited  noise,  sine  and  burst  random.  The  SD-380  communicates 
with  the  computer,  specifically  the  EMODAh  modal  analysis 
software  package,  via  an  IEEE-488  interface,  using  a  General 
Purpose  Interface  Bus  (GPIB) ,  specifically  the  National 
Instruments  AT-GPIB  interface  board.  This  interface  allows 
the  transfer  of  measurement  data  from  the  spectrum  analyzer  to 
the  EMODAL  program,  or  from  the  analyzer  to  an  optional 
digital  pen  plotter. 

In  an  actual  modal  test,  the  SD-380  is  used  to  analyze 
a  vibration  signal,  in  this  case  acceleration,  for  its 
frequency  content.  When  a  structure  is  excited  by  white  noise 
(full  spectrum  excitation),  the  response  is  the  superposition 
of  all  the  modes  of  the  structure.  The  analyzer  determines 
FPEs  from  the  different  measurement  locations  on  the 
structure,  and  performs  a  Fast  Fourier  Transform  to  break  out 
the  separate  modes  at  each  location  on  the  structure.  The 
FRFs  are  shown  as  a  series  of  peaks,  with  characteristic 
center  frequencies,  identifying  them  as  resonances.  It  is 
these  FRFs  that  are  transferred  to  the  EMODAL-PC  software, 
which  creates  a  mathematical  expression  for  each  FRF,  and 
extracts  the  mode  shapes  and  modal  parameters.  The  software 
can  then  calculate  and  display  a  modal  animati.on  for  all 
resonances  of  interest. 
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E.  MODAIi-PC  SOFTWARE  BACRGROUHD 


1 .  Introduction 

with  the  growth  in  popularity  of  modal  analysis  in 
recent  years,  several  good  modal  analysis  software  packages 
have  recently  been  released.  One  of  the  packages  currently  on 
the  market  is  EMODAL-PC  by  Entek  Scientific  Corporation. 
EMODAL-PC  is  a  user-friendly  personal  computer-based  software 
package,  specifically  designed  and  written  for  modal  testing, 
operating  deflection  shape  analysis  and  vibration 
troubleshooting.  It  allows  250  measurement  locations  for  real 
shapes,  125  measurement  locations  for  complex  shapes,  750 
total  DOFs  and  20  mode  shapes. 

EMODAL-PC  requires  an  MS-DOS-based  PC/XT/AT,  80386,  or 
100%  compatible  computer,  with  at  least  a  10  Mb  hard  disk  and 
at  least  640  kb  of  random  access  memory  (of  which  at  least  420 
kb  must  be  available  for  use  by  EMODAL-PC)  .  Other 
requirements  include  at  least  one  floppy  disk  drive,  an  8087, 
80287  or  80387  math  coprocessor,  a  CGA,  EGA  or  VGA  monitor, 
and  an  IEEE-488  interface  port  on  the  computer  on  which  it  is 
to  be  used. 

2 .  Major  Features 

The  version  of  the  EMODAL-PC  software  used  in  the 
Sidewinder  missile  modal  analysis  test  setup  is  version  2.75. 
The  major  features  of  this  version  of  the  EMODAL-PC  software 
include  the  following  [Ref.  3;p.  2-2]: 
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•  Modal  test  outline  to  lead  the  user  through  the  steps 
necessary  to  conduct  a  modal  test. 

•  Easy  editing  of  measurement  location  coordinates  and 
geometry  connectivity. 

•  Automatic  communication  with  a  spectrum  analyzer  via  an 
IEEE-488  interface. 

•  Auto-configuration  of  spectrum  analyzer. 

•  Transfer  of  measured  FRFs  from  the  spectrum  analyzer  to 
the  program. 

•  Three  different  curvefitters ;  co/quad,  circle,  or 
rational  fractional  polynomial  fit. 

•  Real-time  animation  of  deflection  shapes. 

•  Ability  to  plot  undeformed/deformed  geometry. 

•  Ability  to  plot  measurements  in  real,  imaginary, 
magnitude,  phase.  Bode  or  Nyquist  format. 


F.  EMODAL~PC  MODAL  TESTING  FUNDAMENTALS 


The  EMODAL-PC  program  is  a  valuable  tool  for  the 
analyst  or  designer  wishing  to  analyze  and  solve  noise  and 
vibration  problems  of  structures  or  mechanical  systems.  The 
most  successful  approaches  to  solving  these  problems  usually 
utilize  a  combination  of  analysis  and  testing.  In  this  vein, 
the  EMODAL-PC  program  can  be  used  by  the  analyst  in  the 
following  general  sequence  of  steps  [Ref.  3:p.  C-2] . 

•  Step  1:  Identification  of  problem  components  or  systems, 
usually  done  by  operating  vibration  tests. 

•  Step  2:  Perform  experimental  modal  analysis,  in  which 

FRP’s  are  measured  and  the  system's  modal  parameters  are 
identified . 
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•  step  3:  Curvefitting  and  calibration  of  modal  parameters 
to  extract  the  mode  shapes  from  the  frequency  response 
data . 

•  Step  4:  Predictions  of  dynamic  response,  prior  to 

structural  modification  analysis. 

2.  Modal  Testing  vs.  Operating  Deflection  Shape  Test 
a .  Modal  Testing 

Modal  testing  can  be  considered  the  experimental, 
or  measurement-taking  approach,  to  determining  a  structure's 
or  system's  modal  parameters.  In  conducting  a  modal  test,  the 
test  structure  is  subjected  to  some  externally-applied  force. 
An  FRF  is  obtained  (actually  estimated)  from  the  simultaneous 
measurement  of  the  exciting  force  and  the  response.  The  modal 
parameters  are  then  extracted  from  a  set  of  FRITs  obtained  at 
different  locations  on  the  structure. 

Jb.  Operating  Deflection  Shape  T  ist 

An  operating  deflection  shape  test  is  the  process 
of  determining  the  vibration  of  some  structure  or  system  while 
under  operating  conditions.  It  is  not  a  prediction  of  a 
system' s  dynamic  behavior  such  as  would  be  the  result  of  a 
modal  test,  but  instead  it  is  a  direct  measure  of  the  system's 
dynamic  behavior.  To  conduct  an  operating  deflection  shape 
test,  the  structure  (usually  running  machinery)  vibrates  due 
to  self-induced  forces.  The  forces  generated  by  the  structure 
are  not  measured  as  in  modal  testing.  Instead,  the  vibration 
deflection  shape  is  determined  from  a  set  of  vibration 
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amplitude  and  phase  measurements  obtained  at  many  different 
locations  on  the  structure. 

3.  Curvefitting 

In  general  terms,  curvefitting  is  the  process  of 

matching  a  set  of  data  points  obtained  from  FRF  measurements 

I 

.  to  a  mathematical  expression  based  on  a  mass-spring-daunper 

model.  The  si.mplest  type  of  curvefit,  the  straight  line 
curvefit,  is  recognized  by  most  engineers  as  the  result  of  a 
linear  regression  analysis.  In  the  world  of  modal  analysis, 
curvefitting  is  the  first  step  in  the  parameter  estimation 
process  in  which  one  (or  a  few)  measurement  is  used  to 
estimate  the  test  structure's  global  modal  parameters, 

I 

i  specifically  natural  frequency,  damping,  modal  mass  and 

stiffness,  which,  at  least  in  theory,  are  independent  of 
1  position  on  the  structure.  This  is  then  followed  by  the 

second  step  in  the  parameter  estimation  process,  shape 
extraction,  which  uses  all  of  the  measurements  obtained  during 
1  the  test  to  estimate  the  structure's  mode  shapes.  Mode  shape 

I 

coefficients  are  known  as  local  parameters  because  they  do 
depend  on  the  position  on  the  structure. 

a .  Coincident /Quadrature  (Co/Quad)  Fit 

There  are  numerous  curvefitting  techniques 
currently  being  used  in  modal  analysis,  each  having  its 
’  advantages  and  limitations.  The  EMODAL-PC  program  gives  the 

user  a  chioce  of  three  of  the  more  common  curvefitting 
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techniques.  The  simplest  modal  analysis  curvefit  method 
available  is  the  coincident/quadrature  (co/quad)  fit.  This 
method,  which  is  also  called  the  peak  pick  method  because  it 
searches  the  measurement  for  the  spectrum  line  with  the 
maximum  magnitude  in  the  fit  band,  and  uses  this  amplitude  as 
an  estimate  of  the  shape  coefficient,  is  used  mainly  for 
extracting  operating  deflection  shapes,  and  for  extracting 
mode  shapes  quickly  to  check  the  integrity  of  the  data. 
Damping  and  other  modal  parameters  are  not  estimated  in  this 
method.  To  use  this  method,  all  modes  must  be  uncoupled  and 
the  resonance  peaks  must  be  well-separated. 

b.  Circle  Fit 

The  second  type  of  curvefitting  technique  available 
in  EMODAL-PC  is  the  circle  fit.  This  curvefit  method  is  used 
primarily  for  extracting  modal  parameters  when  the  modes  are 
well-separated  and  when  accurate  estimates  of  damping  are  not 
required.  It  performs  a  least  squared  error  algorithm  which 
gives  the  shape  coefficient,  natural  frequency  and  a  rough 
estiiaate  of  the  damping. 

c.  Rational  Fraction  Polynomial  Fit 

The  third  type  of  curvefitting  technique  available 
in  EMODAL-PC  is  the  rational  fraction  polynomial  fit.  This 
curvefit  method  is  a  good,  general  purpose  curvefitter  that 
gives  good  results  in  most  measurement  cases.  Its  main  uses 
include  single  DOF  curvefitting  on  coupled  modes  (situation 


when  two  or  more  modes  are  so  close  in  frequency  that  it.  is 
difficult  to  distinguish  between  peaks)  and  multiple  DOF 
curvefitting.  The  recommended  procedure  to  use  with  this 
curvefit  method  is  to  attempt  to  fit  modes  using  the  single 
DOF  fit,  but  if  the  results  obtained  are  poor,  move  on  to  the 
multiple  DOF  fit.  This  is  the  recommended  curvefitting 
technique  to  use  if  the  most  accurate  results  are  desired. 

G.  EMODAL-PC  TEST  PROCEDURE 

Conducting  a  modal  analysis  of  a  complex  structure  is  not 
necessarily  an  easy  task,  but  the  use  of  modal  analysis 
software  decreases  the  level  of  difficulty.  The  EMODAL-PC 
software  has  a  modal  test  outline  window,  obtained  by  pressing 
<?>  fi'om  the  EMODAL-PC  main  monitor,  which  leads  the  user 
through  the  entire  process  of  setting  up  the  test,  taking 
measurements,  extracting  mode  shapes,  and  .storing  all  the  data 
for  future  use.  The  test  procedure,  as  presented  in  the  modal 
test  outline  [Ref.  13 :p.  2-4],  is  shown  in  Figure  13,  and  is 
defined  in  the  sections  that  follow. 

1 .  Configuring  the  General  Set-up 

The  first  step  in  conducting  a  modal  test  using  the 
EMODAL-PC  software  is  to  confi'-'-  the  software  for  the  test 
about  to  be  conducted.  It  is  in  this  .first  step  that  the  user 
names  the  project  and  tells  the  program  on  what  drive  and 
directory  to  store  the  project  and  measurement  files.  The 
type  of  test  to  be  conducted  (either  modal  test  or  operating 
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Figure  13  EMODAL-PC  Test  Procedure 


deflection  shape  test)  is  chosen,  as  well  as  whether  the 
output  will  be  sent  to  the  computer  screen,  a  printer  or  a 
plotter.  The  most  impoi'tant  input  in  this  step  is  the 
analyzer  GPIB  addi'ess,  without  which  the  program  cannot 
transfer  measurements  from  the  analyzer. 

2 .  Geometry 

The  second  step  of  the  modal  test  outline  is  the  test 
structure  geometry  definition,  in  which  the  user  enters  the 
coordinates  of  the  measurement  locations  on  the  test 
structure,  as  well  as  the  links  (or  elements)  connecting  the 
measurement  locations .  There  are  three  separate  and  necessary 
processes  that  take  place  in  defining  the  structure' s  geometry 
within  the  modal  analy.sis  software.  The  first  is  the  process 
of  entering  the  x,  y  and  z  coordinates  of  all  the  measurement 
locations  into  the  EMODAL-PC  software.  These  measurement 
locations  should  be  completely  thought-out  and  defined  by  the 
analyst  prior  to  entering  them  into  the  computer.  The  add 
function  can  be  used  to  input  additional  measurement 
locations,  should  this  be  necessary,  and  the  edit  function  can 
be  used  to  modify  or  correct  locations  that  have  already  been 
entered.  Once  all  of  the  measurement  locations  are  entered, 
the  next  step  is  to  connect  these  locations  together.  The 
links,  or  elements,  defined  in  this  step  are  the  in.struct ions 
to  the  program  to  accurately  draw  the  test  structure.  Like 
the  coordinates  commands,  the  add  or  insert  function  can  be 
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used  to  er.ter  new  links  and  the  edit  function  is  used  to 


modify  or  coiruect  links  that  have  already  been  entered.  After- 
all  locations  and  links  have  been  entered,  the  display 
geometry  function  should  be  selected  to  verify  that  the 
structure  has  been  drawn  correctly.  This  function  allows  the 
dispJayed  structure  to  be  translated  or  rotated  about  the 
axes,  as  well  as  magnified  to  verify  small,  hidden  or  complex 
part.s  of  the  structure.  Once  the  displayed  geometry  has  been 
verified,  the  analyst  is  ready  to  begin  the  actual  conduct  of 
the  modal  test . 

3 .  Data  Acquisition 

Tile  data  acquisition  process  is  the  most  important 
part  of  a  modal  analysis,  as  only  valrd  measurements  will  give 
valid  test  results.  Obtaining  valid  measurements  requires,  as 
a  minimum,  that  the  analyst  understand  the  operation  of  the 
spectrum  analyzer  in  conjunction  with  the  modal  analysis 
software  package,  that  the  transducers  being  used  have  been 
calibrated  correctly  and  can  measure  response  in  the  frequency 
range  of  interest,  and  that  the  analyst  use  good,  careful  and 
consistent  measurement  technique.  The  correct  procedure  to 
use  in  acquiring  data  is  to  keep  the  reference  transducer 
stationary  throughout  the  test  and  move  the  response 
transducer (s)  to  various  locations  on  the  structure.  The 
spectrum  analyzer  will  then  measure  FRFs  between  the  reference 
and  response  coordinates . 
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The  first  and  most  important  step  in  the  data 
acquisition  process  is  transferring  the  measurements  from  the 
spectrum  analyzer  to  the  program.  To  accomplish  this 
successfully,  the  data  acquisition  parameters  of  the  EMODAL-PC 
software  prog.ram  must  be  configured  properly.  These 
parameters  include  the  project  name,  reference  coordinate, 
measurement  drive  and  directory,  calibration  setup,  and 
measurement  auto-store  feature.  After  the  data  acquisition 
setup  has  been  completed,  the  analyzer  can  be  properly 
configured.  This  can  be  a  rather  complex  task  as  the  SD-380 
has  11  setup  pages,  each  de  licated  to  a  specific  group  of 
functions  or  modes,  and  discussed  in  detail  in  the  SD-380 
Signal  Analyzer  Operator's  Manual  [Ref.  14].  However,  the 
data  acquisition  menu  in  the  modal  test  outline  contains  a 
"Help  setting  up  analyzer  front  panel"  selection,  which 
provides  the  user  with  some  help  in  this  setup  procedure. 
Upon  completion  of  analyzer  setup,  the  "Read  measurement" 
command  can  be  chosen  from  the  data  acquisition  menu,  which 
results  in  the  program  prompting  the  user  for  the  response 
coordinate  for  which  a  transfer  function  is  desired.  After 
the  location  and  direction  (for  example,  lx)  have  been 
entered,  the  program  will  begin  the  data  transfer  process. 
VJhen  completed,  the  response  transducer  (s)  may  be  moved  to  new 
locations,  and  the  process  continued  until  transfer  functions 
have  .en  obtained  for  all  measuremeiit  locations  of  interest. 
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ThG  .oecond  step  in  the  data  acquisition  process  is  the 
storing  and  recalling  of  measurements.  The  recommended  and 
easiest  procedure  for  storing  measurements  is  to  have  the 
program  automatically  store  them  when  they  are  transferred 
from  the  spectrum  analyzer.  This  is  accomplished  via  the 
"Auto-store  measurements=ON"  parameter  in  the  data  acquisition 
setup  menu.  Otherwise,  the  measurements  can  be  stored 
manually  via  the  "Store  measurement"  command  in  the  data 
acquisition  menu.  To  recall  a  measurement  previously  stored, 
simply  select  the  "Recall  measurement"  command  from  the  data 
acquisition  menu,  and  the  program  will  prompt  the  user  for  a 
measurement  number.  The  measurement  number  is  the  numbe.r  the 
program  assigns  to  a  particular  measurement,  and  is  not 
necessarily  the  same  as  the  location  number  on  the  test 
structure  at  which  the  measurement  was  made.  To  determine  the 
measurement  number  assigned  to  a  particular  measurement, 
select  the  "Measurement  directory"  command  from  the  data 
acquisition  menu,  enter  the  project  name  at  the  prompt,  then 
enter  "a"  to  obtain  a  directory  listing  of  all  measurements. 
When  the  m.easurement  number  of  interest  is  determined,  enter 
that  nunu^er  at  the  "Recall  measurement"  prompt,  and  the 
program  will  then  recall  that  measurement. 

The  third  step  in  the  data  acqiiisition  process  is 
viewing  measurement,  usually  performed  after  recalling  a 
measurement.  lifter  selecting  the  "View  measurement"  command 
from  the  data  acquisition  menu,  the  program  displays  the 
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recalled  measurement,  along  with  commands  which  allow  the  user 
to  magnify  or  expand  the  view,  set  up  frequency  bands  for 
shape  extraction  and  execute  curvefits . 

4 .  Extract  and  Animate  Shapes 

This  is  the  step  in  a  modal  test  where  fit  bands  are 
set  up  and  curvefits  are  perforir.ed,  and  mode  shapes  are 
extracted  and  animated.  The  purpose  of  setting  up  fit  bands 
is  to  identify  the  frequencies  of  interest  on  a  displayed 
measurement  for  which  a  curvefit  wil.l  be  conducted.  The 
procedure  for  setting  up  fit  bands  and  conducting  a  curwefit 
is  simple  and  ctraight-forward,  with  the  program  prompting  the 
use."  for  the  desired  input.  From  the  "View  measurement  and 
set  up  fit  bands"  option  in  the  "Extract  and  Animate  Shapes" 
menu,  the  <SETUP  FIT>  option  should  be  selected.  The  program 
will  prompt  the  user  for  a  shape  number  (assigned  in  order  of 
increasing  frequency) ,  curvefit  type  (either  co/quad,  circle 
or  polynomial)  and  the  upper  and  lower  limits  of  the  desired 
curvefit  band.  The  user  then  selects  the  <CURVEFIT>  option  to 
execute  the  curvefit,  for  which  the  program  displays  a  graph 
of  the  curvefit  results  in  real,  imaginary,  log  magnitude  and 
Nyquist  formats.  After  the  curvefit,  the  user  has  the  option 
of  rejecting  the  curvefit  or  accepting  it,  if  the  results  are 
satisfactory,  and  saving  the  results  in  the  project's 
parameter  table . 


61 


To  extract  shapes,  each  measurement  on  the  disk  is 
curvefit  using  the  above  process.  For  each  fit,  a  shape 
coefficient  is  computed  and  sa'^’ed  in  the  project's  shape  file, 
from  which  the  structure's  mode  shapes  are  generated.  To 
initiate  the  shape  extraction  process,  select  the  "Extract 
shapes"  option  from  the  "Extract  and  Animate  Shapes"  menu, 
then  select  tne  "Shape  Extraction  Setup"  option  and  set  the 
parar-eters  as  desired.  The  parameters  include  extract  type 
(real  or  complex) ,  reference  coordinate,  and  measurement 
drive/directory.  After  these  parameters  have  been  set 
properly,  select  the  "Extract  shape.s"  option,  and  enter  the 
measurement  numbers  for  which  extracted  shape  coefficients  are 
desired.  The  program  will  recall  these  measurements  from  disk 
and  curvefit  them,  resulting  in  the  generation  of  the 
structure's  mode  shapes,  which  can  then  be  animated. 

The  purpose  of  animating  the  shapes  is  to  observe  how 
the  different  measurement  locations  on  the  structure  move, 
either  in  translation  or  rotation  or  both,  as  a  result  of  the 
excitation  force.  After  the  shapes  have  been  extracted,  the 
user  selects  tne  "Animate  shapes"  option  from  the  "Extract  and 
Animate  Shapes"  menu,  and  enters  the  desired  shape  number. 
The  program  then  displays  tne  animation  of  the  deformed 
structure.  Numerous  options  are  given  in  the  animation  menu, 
including  the  ability  to  zoom  into  different  parts  of  the 
structure,  rotate  or  tx-anslate  the  geometry,  and  display  both 
the  deformed  and  undeformed  geometry  on  dual  plots. 
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5.  Store  Project 

All  of  the  files  associated  with  a  particular  project 
can  be  stored  very  easily  by  selecting  the  "Store  Project" 
menu  from  the  modal  test  outline.  It  is  very  important  to 
store  the  project  files  immediately  upon  completion  as  the 
files  will  be  lost  if  the  program  stops  for  any  reason  (e.g., 
power  outage  or  user  exiting  the  program  without  saving) . 
There  are  actually  six  different  types  of  project  files  that 
are  stored  when  this  option  is  selected.  Each  of  these  files, 
when  saved,  is  given  the  same  name  as  the  project  name  chosen 
when  the  general  setup  is  configured,  but  each  has  a  different 
filename  extension,  as  summarized  below  [Ref.  13;p.  2-61]. 

•  Coordinates  file  -  contains  the  x,  y,  and  z  coordinates  of 
each  measurement  location;  has  extension  .CRD 

•  Lii\ks  file  -  contains  the  connectivity  that  draws  the  test 
structure;  has  extension  .LNK 

»  Parameter  table  file  -  contains  the  fit  band  definitions 
and  the  modal  parameters  extracted  by  curvefitting;  has 
extension  .PAR 

•  Shapes  file  -  contains  the  deflection  shapes,  arranged  in 
columns  by  shape  number;  has  extension  .Sxx 

»  Setup  menus  file  -  contains  program  configuration 
parameter  settings;  las  extension  .SET 

•  Title  file  -  contains  uhe  project  title;  has  extension 
.  TTL 

Once  all  of  the  project  files  have  been  stored  on  the 
computer's  hard  disk,  they  can  easily  be  copied  on  floppy 
disks  for  future  reference. 


6.  Recall  Project 

A  previously  stored  project  may  need  to  be  recalled 
back,  into  computer  memory  for  the  purpose  of  reviewing  the 
results  or  preparing  reports.  To  recall  a  project,  simply 
select  the  "Recall  Project"  step  of  the  modal  test  outline  and 
enter  the  project  name  at  the  prompt.  The  program  will  recall 
all  of  the  project  files  and  upon  completion,  will  redisplay 
the  modal  test  outline. 

7 .  Prepare  Reports 

This  final  step  allows  the  user  to  print  project 
files,  produce  plots  of  the  undeformed  and  deformed  structure, 
and  produce  plots  of  specific  measurements,  all  in  report- 
ready  format.  Prior  to  receiving  any  of  these  outputs  from 
the  program,  however,  the  program  must  be  correctly  configured 
for  the  printer  and/or  plotter  to  be  used.  This  is 
accomplished  during  the  "General  Program  Setup"  step.  After 
verification  of  the  correct  configuration,  and  upon  selection 
of  the  "Prepare  Reports"  option  of  the  modal  test  outline,  the 
user  is  given  four  options.  The  first  option,  "Print  all 
files  and  tables",  displays  the  print  selection  menu,  from 
which  the  file  or  files  ho  be  printed  are  selected.  The 
second  option,  "Plot  shapes",  will  produce  plots  of  either  the 
undeformed  or  deformed  structure,  sending  these  plots  either 
to  the  computer  screen,  or  to  an  external  printer  or  plotter. 
The  third  option,  "Recall  a  measurement  from  disk",  is  used  to 
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recall  a  specific  measure- \ent  from  the  measurement  directory 
for  the  purpose  of  plotting  it.  The  fourth  option,  "Plot 
measurement",  can  then  be  selected  to  produce  a  plot  of  that 
measurement .  However,  the  "Measurement  plotting  setup"  menu 
must  first  be  configured  for  plotting  the  measurement  of 
interest.  In  this  menu,  the  user  selects  plotting  parameters 
including  type  of  plot  (real,  imaginary,  Bode,  etc.),  line 
type,  title,  x— axis  and  y-axis  labels  and  divisions,  scaling, 
and  so  on. 
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IV.  CONCLUSIONS 


A.  ANALYTICAL  MODELING 

Within  the  scope  of  this  investigation,  the  MSC/PAL2 
finite  element  software  is  satisfactory  for  obtaining  an 
analytical  model  of  a  missile  system,  and  analytically 
determining  its  natural  frequencies  and  mode  shapes .  The 
conduct  of  a  normal  modes  analysis  yields  insight  into  the 
dynamic  behavior  of  the  structure  under  a  wide  variety  of 
dynamic  loading.  Knowing  the  lowest  few  resonant  frec[uencies 
and  mode  shapes  provides  the  designer  or  analyst  with 
information  which  may  aid  in  identifying  potential  areas 
requiring  redesign  or  possible  interface  problems  between 
components  of  a  system. 

In  the  case  of  the  AIM-9  missile,  the  material  properties 
and  mass  distributions  in  the  model  obtained  via  the  PAL2 
software  can  be  modified,  based  on  the  outcome  of  an 
experimental  comprehensive  modal  test,  to  accurately  represent 
the  dynamic  response  of  the  actual  missile.  Once  an  accurate 
model  is  obtained,  feasibility  studies  can  easily  be  performed 
for  the  AIM-9  upgrades  curreni^ly  being  planned.  Specifically, 
the  AIM-9M-8/9  program  will  retrofit  existing  missiles  with  an 
infrared  counter-countermeasures  feature,  the  AIM-9R  will 
incorporate  an  imaging  focal  plane  array  to  increase  target 
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acquisition  range  while  improving  the  missile's  counteirmeasure 
and  clutter  rejection  capabilities,  and  the  year  2000  AIM-SX 
will  feature  an  improved  seeker,  airframe,  and  warhead/fuze 
[Ref.  15]  .  All  of  these  changes  may  affect  the  material's 
parameters  and  mass  distributions,  which  can  be  readily 
modified  in  the  finite  eleme..it  model. 

The  AIM-9  model  developed  in  this  investigation 
satisfactorily  represents  an  actual  missile;  however,  the 
material's  parameters  and  mass  distribution  have  yet  to  be 
matched  to  the  natural  frequencies  and  mode  shapes  deuermined 
by  the  dynamic  modal  measurement  of  an  actual  missile.  Modal 
tests  are  planned  for  a  follow-on  phase  of  the  investigation, 
at  which  point  the  model  can  eas.rly  be  modified  to  reflect  an 
accurate  dynamic  response  by  altering  the  mass  locations, 
element  stiffness  or  material  properties. 

B .  EXPERIMENTAL  MODELING 

While  the  GVT  stand  is  satisf  act  ..^ry  for  longitudinal 
excitation  of  a  test  article,  it  is  presently  unsuitable  for 
excitation  of  a  test  article  in  the  lateral  direction  and  must 
be  structurally  stiffened  prior  to  the  conduct  of  a 
comprehensive  mi.ssile  modal  test. 

Within  the  scope  of  this  investigation,  the  modal  analysis 
test  setup  and  the  associated  methodologies  established  are 
satisfactory  for  conducting  a  comprehensive  modal  analysis  of 
th  AIM-9  missile  system  for  the  purpose  of  experimentally 
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determining  the  missile's  transient  response  and  modal 
parameters.  The  Fast  Fourier  Transform-based  analyzer  and  the 
EMODAL-PC  modal  analysis  software  have  been  incorporated  into 
the  test  setup,  and  are  ready  for  transfer  function 
measurement  and  subsequent  determination  of  a  structure's 
modal  parameters . 

C.  COMBINED  TESTING  AND  ANALYSIS 

The  use  of  the  analytical  and  experimental  techniques  in 
a  complementary  manner  allows  the  designer  or  analyst  to  draw 
on  the  strengths  of  each  technique  to  thoroughly  and 
conpletely  describe  the  dynamic  response  of  a  structure.  By 
first  determining  a  structure' s  dynamic  response 
experimentally  and  then  modifying  the  parameters  in  an 
analytical  model  to  match  that  response,  the  analyst  can  use 
this  analytical  model  to  investigate  different  design  options, 
sometimes  called  "what  if?"  studies,  thereby  avoiding  the 
costly  and  time-consuming  trial— and— error  approach  to  solving 
vibration  problems  on  the  actual  hardware.  The  "what  if?" 
approach  can  be  iterated  on  the  model  before  any  hardware 
modifications  or  upgrades  are  implemented. 
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V.  RECOMMENDATIONS 


A.  ANALYTICAL  MODELING 

After  completion  of  a  modal  test  on  the  AIM-9  missile 
system,  the  analytical  model  should  be  modified  to  match  the 
same  mode  shapes  and  natural  frequencies  observed 
expeirimentally .  This  can  be  accomplished  by  changing  the 
missile's  materral  properties,  by  re-distributing  the  lumped 
masses  more  evenly  along  the  missile  body,  by  adding  stiffness 
to  the  front  fins  and  by  defining  the  missile' s  attachment 
points  across  an  element  vice  at  a  sing].e  node  point.  The 
final  and  refined  version  of  the  analytical  model,  and  hence 
the  actual  modifications  made  to  the  model,  v;ill  depend 
entirely  on  the  mode  shapes  and  frequencies  obtained  via  modal 
analysis . 

Upon  completion  of  an  accurate  AIM-9  missile  analytical 
model,  design  iterations  should  be  conducted  for  the  planned 
AIM-9  upgrades  to  analytically  determine  the  natural 
frequencies  and  mode  shapes  of  these  upgrades,  in  an  effort  to 
ascertain  whether  a  flutter  problem  will  exist  between  the  P-3 
wing  (or  the  Orion  II  wing)  and  the  upgraded  AIM-9  missile. 
This  can  be  accomplished  primarily  by  modifying  the  model's 
mass  distribution  for  the  AIM-9M-8,/9  and  AIM-9R,  and  by 
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modifying  the  model's  mass  distribution,  fin  surfaces  and 
possibly  the  shape  of  the  body  for  the  AIM-9X  version. 

B .  EXPERIMENTAL  MODELING 

The  GVT  stand  needs  to  be  stiffened  laterally  to  allow 
excitation  measurements  in  the  lateral  direction.  The 
structure  could  be  stiffened  in  several  ways,  including 
replacing  the  upper  and  lower  3/16"  plates  by  thicker  1/2" 
plates,  adding  sandbags  between  the  plates  to  provide  damping, 
and  adding  transverse  bracing  between  the  three  internal  "I" 
beams . 

The  AIM-9  missile  system  should  undergo  a  comprehensive 
modal  test  to  completely  determine  its  mode  shapes  and  natural 
frequencies  experimentally  in  both  the  longitudinal  and 
lateral  directions,  using  the  missile  modal  analysis  test 
setup  completed  in  this  investigation.  The  results  of  this 
comprehensive  modal  test  should  then  be  compared  to  the  P-3 
Orion  and  Orion  II  wing  models  to  .examine  any  potential 
flutter  problems  at  any  of  the  proposed  wing  carriage 
stations . 

C.  COMBINED  TESTING  AND  ANALYSIS 

The  facilities  and  methodologies  discussed  and  used  in 
this  investigation  of  the  AIM— 9  missile's  dynamic  response 
should  be  utilized  to  conduct  similar  combined  testing  and 
analysis  on  other  candidate  missile  systems,  including  the 
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AGM-65  Maverick  and  AGM-88  HAIJvI,  considered  for  incorporation 
into  the  P-3's  weapon  inventory. 
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APPENDIX  A 


The  MSC/PAL2  software  code  for  the  finite  element  model  of 
the  AIM— 9  missile  is  given  below.  Explanatory  comment 
statements  W'ere  included  where  appropriate.  For  those  PAL2 
commands  which  could  not  be  fully  explained  in  a  brief  comment 
statement,  the  appropriate  page  of  the  MSC/PA’j2  User '  s  Manual 
[Ref.  6]  is  given,  on  which  a  detailed  discussion  of  the 
command  can  be  found. 


TITLE  -  AIM9  GEOMETRY  FILE 
C 

C  DEFINE  NODES  FOR  MISSILE  BODY  IN  CYLINDRICAL  COORDINATES 
C 

NODE  .3 

C  NODE  NUMBER,  RADIUS,  THETA,  Z 
C  (Ref.  b:p.  3-92) 

1  2.50,  0,  0  THROUGH  16  2.50,  337.5,  0 

545  2.50,  0,  102  THROUGH  560  2.50,  337.5,  102 

561  0,  0,  0 

C 

C  GENERATE  NODES  OF  MISSILE  BODY  AUTOMATICALLY 


C  (Ref, 

NODE  22 

,  6 

:p.  3-99) 

1 ,  16,  .560, 

1, 

16 

NODE  22 

1 ,  ] 6 ,  560, 

15, 

.  16 

C 

C  GENERATE  NODES  OF  TAPERED  NO5E0ONE  AUTOMATICALLY 
C 

NODE  3 

562  2.50,  0,  103  THROUGH  577  2.50,  337.5,  103 
594  1.25,  0,  112  THROUGH  609  1.25,  337.5,  112 

NODE  22 

562,  577,  bC9,  1,  16 
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NODE  22  ,  . 

562,  577,  609,  15,  16 

NODE  3 

610  0,  0,  113 
C 

C  CHANGE  POSITION  OF  CERTAIN  NODES  TO  ALLOW  ATTACHMENT  OF  FINS 
C 

C  MOVE  SECOND  S  TENTH  CONCENTRIC  RINGS  1  INCH  IN  -Z  DIRECTION 

C 

NODE  1 

17  0,  0,  -1  THROUGH  32  0,  0,  -1 

NODE  1 

145  0,  0,  -1,  THROUGH  160  0,  0,  -1 

C  MOVE  34TH  CONCENTRIC  RING  2  INCHES  IN  THE  -Z  DIRECTION 
NODE  1 

529  0,  0,  -2,  THROUGH  544  0,  0,  -2 
C 

C  DEFINE  NODE  POINTS  FOR  FRONT  FINS 
C 

NODE  3 
C 

C  FIRST  FRONT  FIN 
C 


611 

2.75, 

45,  93 

612 

2.75, 

45,  96 

613 

2.75, 

45,  97 

614 

2 . 75, 

45,  101.5 

615 

7.00, 

45,  97.25 

616 

7 . 00, 

45,  95 

617 

7 . 00, 

45,  94.5 

618 

7.00, 

45,  93 

619 

11.25 

,  45,  93 

C 

C  SECOND 

FRONT  FIN 

c 

620 

2.75, 

135,  93 

621 

2.75, 

135,  96 

622 

2.75, 

135,  97 

623 

2 . 75, 

135,  101.5 

624 

7.00, 

135,  97.25 

625 

7.00, 

135,  95 

626 

7.00, 

135,  94.5 

627 

7.00, 

135,  93 

628 

11 . 25 

,  135,  93 
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c 


C  THIRD  FRONT 

C 

FIN 

629 

2.75, 

225, 

93 

630 

2.75, 

225, 

96 

631 

2,75, 

225, 

97 

632 

2.75, 

225, 

101 .5 

633 

7.00, 

225, 

97 . 25 

634 

7  .  CO, 

225, 

95 

635 

7.00, 

22  5, 

94.5 

636 

7.00, 

225, 

93 

637 

C 

11.25 

,  225 

,  95 

C  FOURTH 

C 

FRONT 

FIN 

638 

2.75, 

315, 

93 

639 

2.75, 

315, 

96 

640 

2.75, 

315, 

97 

641 

2.75, 

315, 

101.5 

642 

7.00, 

315, 

97 . 25 

643 

7.00, 

315, 

95 

644 

7.00, 

315, 

94 . 5 

645 

7.00, 

315, 

93 

646 

11.25 

,  315 

,  93 

c 

C  DEB'INE  NODE  POINTS  FOR  REAR  WINGS 
C 

NODE  3 
C 


C  FIRST  REAR  WING 
C 


64  7 

4  . 

75, 

45, 

2 

648 

4  . 

75, 

45, 

6 

649 

4  . 

75, 

45, 

9 

650 

4  . 

75, 

45, 

12 

651 

1  . 

75, 

45, 

15 

652 

4  . 

75, 

45, 

18 

653 

4  . 

75, 

45, 

21 

654 

4  . 

75, 

45, 

24 

655 

7  . 

88, 

45, 

21 

656 

11 

.  50, 

45, 

18 

657 

li 

.50, 

45 

15 

658 

11 

.50, 

45, 

12 

659 

11 

.  50, 

45, 

9 

660 

11 

.50, 

45, 

6 

661 

11 

.50, 

45, 

2 

C 

C  SECOND  RF.AR  WING 
C 

662  4.75,  135,  2 
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I 


663 

4  . 

75, 

135, 

6 

6  64 

4  . 

75, 

135, 

9 

665 

4  . 

75, 

135, 

12 

666 

4  . 

75, 

135, 

15 

667 

4  . 

75, 

135, 

18 

668 

4  . 

75, 

135, 

21 

669 

4  . 

75, 

135, 

24 

670 

7  . 

88, 

135, 

21 

671 

11 

.50, 

135, 

.  18 

672 

11 

.50, 

135, 

15 

673 

11 

.50, 

135, 

12 

674 

11 

.50, 

135, 

9 

675 

11 

.50, 

135, 

6 

676 

11 

.50, 

135, 

2 

c 


C  THIRD  REAR  WING 
C 


67  7 

4.75, 

225,  2 

678 

4.75, 

225,  6 

679 

4.75, 

225,  9 

680 

4.75, 

225,  12 

681 

4.75, 

225,  15 

682 

4.75, 

225,  18 

683 

4.75, 

225,  21 

684 

4.75, 

225,  24 

685 

7.88, 

225,  21 

68  6 

11 .50, 

0  O  C  IQ 

^  ^  ^  f  -1.  V 

687 

11.50, 

225,  15 

688 

1]  .50, 

225,  12 

689 

11.50, 

225,  9 

690 

11.50, 

225,  6 

691 

11.50, 

225,  2 

C 

C  FOURTH  REAR  WING 

C 

692 

4 .75, 

315,  2 

693 

4.75, 

315,  6 

694 

4.75, 

315,  9 

695 

4 . 75, 

315,  12 

696 

4 . 75, 

315,  15 

6  97 

4 . 75, 

315,  18 

698 

4.75, 

315,  21 

699 

4 . 75, 

315,  24 

700 

7.88, 

315,  21 

701 

11.50, 

315,  18 

702 

11.50, 

315,  15 

703 

11.50, 

315,  12 

704 

11 . 50, 

315,  9 

705 

11.50, 

315,  6 

706 

11 . 50, 

315,  2 

75 


c 

C  ENTER  r^ATERIAL  PROPEP.TIES 
C 

C  YOUNG'S  MODULUS,  SHEAR  MODULUS,  I-IASS  DENSITY,  POISSON'S 
C  RATIO,  TENSILE  YIELD  STRESS 
C  (Ref.  6:p.  3-83) 

MAT  30.0E6,  12,0E6,  7.76E-4,  0.25,  30 . 0E3 
C 

C  DEFINE  PLATE  TYPE  FOR  MISSILE  BODY 
C 

C  QUADRILATERAL  PLATE  ELEMENT,  Nl,  N2 ,  THICKNESS 
C  (Ref.  C:p.  3-113) 

QUAD  1,  1,  0.25 
C 

C  GENERATE  QUAD  PLATE  ELEMENTS  ON  MISSILE  BODY  AUTOMATICALLY 
C  (Ref.  6:p.  3-.56) 

GENERATE  CONNECTS  1,  16,  560,  1,  16 

GENERATE  CONNECTS  1,  16,  560,  15,  16 
GENERATE  CONNECTS  562,  577,  609,  1,  16 

GENERATE  CONNECTS  562,  577,  609,  15,  16 

C 

C  DEFINE  REAR  MISSILE  PLATE  AND  CONNECTIVITY 
C 

C  TRIANGULAR  PLATE  ELEMENT,  Nl,  N2,  THICKNESS 
C  (Ref.  6:p.  3-121) 

TRI  1,  2,  0.25 

Q 

C  CONNECT  NODE  A  TO  NODE  B  TO  NODE  C 
C  (Ref.  6;p.  3-27) 

CON  1  TO  2  TO  561 
CON  2  TO  3  TO  561 
CON  3  TO  4  TO  561 
CON  4  TO  5  TO  561 
CON  5  TO  6  TO  561 
CON  6  TO  7  TO  561 
CON  7  TO  8  TO  561 
CON  8  TO  9  TO  561 
CON  9  TO  10  TO  561 
CON  10  TO  11  TO  561 
CON  11  TO  12  TO  561 
CON  12  TO  13  TO  561 
CON  13  TO  14  TO  561 
CON  14  TO  15  TO  561 
CON  15  TO  16  TO  561 
CON  16  TO  1  TO  561 
C 

C  DEFINE  PLATE  CAP  CONNECTIVITY  OF  NOSECONE 
C 


CON 

594 

TO 

5  95 

TO 

610 

CON 

595 

TO 

596 

TO 

610 

CON 

596 

TO 

597 

TO 

610 
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CON 

597 

TO 

598 

TO 

610 

CON 

598 

TO 

599 

TO 

610 

CON 

599 

TO 

60C 

TO 

610 

CON 

600 

TO 

601 

TO 

610 

CON 

603 

TO 

602 

TO 

610 

CON 

602 

TO 

603 

TO 

610 

CON 

603 

TO 

604 

TO 

610 

CON 

604 

TO 

605 

TO 

610 

CON 

605 

TO 

606 

TO 

610 

CON 

606 

TO 

607 

TO 

610 

CON 

60  7 

TO 

608 

TO 

610 

CON 

608 

TO 

609 

TO 

610 

CON 

609 

TO 

594 

TO 

610 

c 

C  CONNECT  BODY  TO  NOSSCONE 
C 

QUAD  1,  1,  0.25 
C 


C  CONNECT  NODE  A  TO  NODE  E  TO  NODE  C  TO  NODE  D 
C  (CONNECT  COUNTERCLOCKWISE) 


C 

(Ref . 

6;p. 

,  3-2 

7) 

CON 

545 

TO 

562 

TO 

563 

TO 

546 

CON 

546 

TO 

563 

TO 

564 

TO 

547 

CON 

547 

TO 

564 

TO 

565 

TO 

548 

CON 

548 

TO 

565 

TO 

566 

TO 

549 

CON 

549 

TO 

566 

TO 

567 

TO 

550 

CON 

550 

TO 

567 

TO 

568 

TO 

551 

CON 

551 

TO 

568 

TO 

569 

TO 

552 

CON 

552 

TO 

569 

TO 

570 

TO 

553 

CON 

553 

TO 

570 

TO 

571 

TO 

554 

CON 

554 

TO 

571 

TO 

572 

1  ■) 

555 

CON 

555 

TO 

572 

TO 

573 

T  > 

556 

CON 

556 

TO 

573 

TO 

574 

TO 

557 

CON 

557 

TO 

574 

TO 

575 

TO 

558 

CON 

558 

TO 

575 

TO 

576 

TO 

559 

CON 

559 

TO 

576 

TO 

577 

TO 

560 

CON 

560 

TO 

577 

TO 

562 

TO 

545 

C 

C  DEFINE  FRONT  FIN  PLATE  TYPE  AND  CONNECTIVITY 
C 

C  FIRST  FIN 
C 

C  QUAD  PLATE  TO  ATTACH  FIN  TO  MISSILE  BODY 
C 

QUAD  1,  1,  0.35 

CON  515  TO  531  TO  613  TO  612 

C 

C  FIN  CONSTRUCTION 

C 

QUAD  1 ,  1 ,  0.10 

CON  611  TO  612  TO  617  TO  618 
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CON 

612 

TO 

613 

TO 

616 

TO 

617 

CON 

613 

TO 

614 

TO 

615 

TO 

616 

TRI 

1,  2 

,  0 

.10 

CON 

619 

TO 

61  8 

TO 

617 

CON 

619 

TO 

617 

TO 

616 

CON 

619 

TO 

616 

TO 

615 

C 

C  SECOND 

i  FRONT 

FIN 

C 

QUAD 

1, 

1, 

0 . 35 

CON 

519 

TO 

535 

TO 

622 

TO 

621 

QUAD 

1, 

1, 

0 . 1 

CON 

620 

TO 

621 

TO 

626 

TO 

627 

CON 

621 

TO 

622 

TO 

625 

TO 

626 

CON 

622 

TO 

623 

TO 

624 

TO 

625 

TRI 

1,  2 

,  0 

.  1 

CON 

628 

TO 

627 

TO 

626 

CON 

628 

TO 

626 

TO 

625 

CON 

628 

TO 

625 

TO 

624 

C 

C  THIRD 

FRONT  FIN 

C 

QUAD 

1, 

1, 

0 . 35 

CON 

523 

TO 

539 

TO 

631 

TO 

630 

QUAD 

1, 

1, 

0 . 1 

CON 

629 

TO 

630 

TO 

635 

TO 

636 

CON 

630 

TO 

631 

TO 

634 

TO 

635 

CON 

631 

TO 

632 

TO 

633 

TO 

634 

TRI 

1,  2 

,  0 

.  1 

CON 

637 

TO 

636 

TO 

635 

CON 

637 

TO 

635 

TO 

634 

CON 

637 

TO 

634 

TO 

633 

C 

C  FOURTH 

:  FRONT 

FIN 

C 

QUAD 

1, 

1, 

0 . 35 

CON 

52  7 

TO 

543 

TO 

640 

TO 

639 

QUAD 

1, 

1, 

0 . 1 

CON 

638 

TO 

639 

TO 

644 

TO 

645 

CON 

639 

TO 

640 

TO 

643 

TO 

644 

CON 

640 

TO 

641 

TO 

642 

TO 

643 

TRI 

1,  2 

,  0 

.  JL 

CON 

646 

TO 

645 

TO 

644 

CON 

646 

TO 

64  4 

TO 

643 

CON 

645 

TO 

643 

TO 

642 

c 

C  DEFINE  iyiATERIA.L  PROFEP.TIES  FOR  REAR  WINGS 
C 

MAT  10E6,  4E6,  2.59E-4,  0.25,  10E3 
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c 

C  DEFINE  REAR  WING  PLATE  TYPE  AND  CONNECTIVITY 
C 

C  FIRST  REAR  WING 
C 

QUAD  1 ,  1 ,  0.5 
CON  19  TO  35  TO  648  TO  647 
CON  647  TO  648  TO  660  TO  661 
CON  35  TO  51  TO  649  TO  648 
CON  648  TO  649  TO  659  TO  660 
CON  51  TO  67  TO  650  TO  649 
CON  649  TO  650  TO  658  TO  659 
CON  67  TO  83  TO  651  TO  650 
CON  650  TO  651  TO  657  TO  658 
CON  83  TO  99  TO  652  TO  65: 

CON  651  TO  652  TO  656  TO  657 
CON  99  TO  115  TO  653  TO  652 
CON  652  TO  653  TO  655  TO  656 
CON  115  TO  131  TO  654  TO  653 
TRI  1,  2,  0.5 
CON  654  TO  131  TO  147 
CON  655  TO  653  TO  654 
C 

C  SECOND  REAR  WING 
C 

QUAD  1 ,  1 ,  0.5 

CON  23  TO  39  TO  663  TO  662 


CON 

6  62 

TO  663 

TO 

675 

TO 

67  6 

CON 

39 

TO  55  TO 

i  664  TO 

'  663 

CON 

663 

TO  664 

TO 

674 

TO 

675 

CON 

55 

TO  71  TO 

i  665  TO 

'  664 

CON 

664 

TO  665 

TO 

673 

TO 

674 

CON 

71 

TO  87  TO 

i  666  TO 

'  665 

CON 

665 

TO  666 

TO 

672 

TO 

673 

CON 

87 

TO  103  TO  6 

67  TO  ' 

666 

CON 

666 

TO  667 

TO 

671 

TO 

67^^ 

CON 

103 

TO  119 

TO 

668 

TO 

66  ; 

CON 

667 

TO  668 

TO 

670 

TO 

671 

CON 

119 

TO  135 

TO 

669 

TO 

668 

TRI 

1, 

2,  0.5 

CON 

669 

TO  135 

TO 

151 

CON 

670 

TO  668 

TO 

669 

C 

C  THIRD  REAR  WING 
C 

QUAD  1 ,  1 ,  0.5 
CON  27  TO  43  TO  678  TO  677 
CON  677  TO  676  TO  690  TO  691 
CON  43  TO  59  TO  679  TO  678 
CON  678  TO  679  TO  689  TO  690 
CON  59  TO  75  TO  680  TO  679 
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CON 

679 

TO 

680 

TO 

688 

TO 

689 

CON 

75 

TO 

91  TO  68 

m 

j. 

0  680 

CON 

680 

TO 

681 

TO 

V.1 8  7 

TO 

688 

CON 

91 

TO 

107 

TO  682 

TO 

681 

CON 

681 

TO 

682 

TO 

686 

TO 

687 

CON 

107 

TO 

123 

TO 

683 

TO 

682 

CON 

682 

TO 

683 

TO 

685 

TO 

686 

CON 

123 

TO 

139 

TO 

684 

TO 

683 

TRI 

1, 

2, 

0 . 5 

CON 

684 

TO 

139 

TO 

155 

CON 

685 

TO 

683 

TO 

684 

c 

C  FOURTH  REAR  WING 
C 

QUAD  1 ,  1 ,  0.5 


CON 

31 

TO  47  TO 

693  TO 

692 

CON 

692 

TO  693 

TO 

705 

TO 

706 

CON 

47 

TO  63  TO 

694  TO 

693 

CON 

693 

TO  694 

TO 

704 

TO 

705 

CON 

63 

TO  79  TO 

6  95  TO 

694 

CON 

694 

TO  695 

TO 

703 

TO 

704 

CON 

79 

TO  95  TO 

696  TO 

695 

CON 

695 

TO  696 

TO 

702 

TO 

703 

CON 

95 

TO  111  TO 

697  TO 

696 

CON 

696 

TO  697 

TO 

701 

TO 

702 

CON 

111 

TO  127 

TO 

698 

TO 

697 

CON 

697 

TO  698 

TO 

700 

TO 

701 

CON 

/ 

TO  14  3 

TO 

699 

TO 

r» 

o  ::?o 

TRI 

1, 

2,  0.5 

CON 

699 

TO  143 

TO 

159 

CON 

700 

TO  698 

TO 

699 

C 

C  DEFINE  ATTACHMENT  POINTS 
C 

C  AT  ATTACHMENT  POINTS,  ALLOW  ROTATION  ABOUT  X  AXIS  ONLY 


C 

(Ref 

.  6 

ZERO  1 

TA 

85, 

213, 

373 

RY 

85, 

213, 

373 

RZ 

35, 

213, 

373 

C 

C  DEFINE  POINT  MASSES  FOR  SEPARATE  MISSILE  COMPONENTS 
C 

C  MASS  COMMAND,  NODE  NUMBER,  MASS  OF  COMPONENT 
C 

C  MASS  OF  MOTOR  COMPONENT  (TWO  POINT  MASSES  SINCE  70"  LONG) 
C 

MASS  141,  0.128 
MASS  269,  0.128 


80 


APPENDIX  B 


Figures  14  through  19  depict  the  first  six  mode  shapes  for 
the  AIM-9  missile  model  constructed  in  this  investigation. 
They  represent  the  missile's  deformations  at  its  first  six 
natural  frequencies  of  29.65,  50.79,  52.97,  57.36,  57.89,  and 
58.80  hz .  It  should  be  noted  that  these  frequencies,  and 
hence  these  mode  shapes,  are  not  necessarily  representative  of 
the  actual  missile  response.  These  results  were  obtained 
mathematically,  strictly  from  the  model  geometry,  mass 
density,  and  additi.onal  lumped  masses  representing  the 
individual  missile  components.  The  figures,  presented  in  the 
optimum  planar  representation  for  viewing  the  reponses,  show 
still  frames  of  the  animated  mode  shapes;  the  first  view 
showing  the  upper  limit  of  deformation,  the  second  showing  the 
undeformed  structure,  and  the  third  showing  the  lower  limit  of 
def ormatioii .  While  not  totally  obvious  from  the  figures,  the 
missile  model  exhibited  the  following  modal  responses. 

•  Mode  1  (29.65  hz)  -  whole  body  z-axis  rotation  about 

attachment  points  with  some  bending  in  x-z  plane 

•  Mode  2  (50.79  hz)  -  bending  mode  in  y-z  plane  with  slight 
front  fin  z-axis  i'otation 

•  Mode  3  (52.97  hz)  -  bending  mode  in  x-z  plane  with  slight 
front  fin  z-axis  rotation 

•  Mode  4  (57.36  hz)  -  front  fin  z-axis  rotation  only;  no 

body  motion 

•  Mode  5  (i  ■  89  hz)  -  front  fin  z-axis  rotation  with  some 

body  torsion  in  the  x— y  plane 
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•  Mode  6  (58.80  hr;)  -  front  fin  r— axis  rotation  with  some 

body  bending  in  the  y-z  plane  and  torsion  in  the  x-y  plane 

The  figures,  showing  what  app'ear  to  be  strange  contortions 

of  the  front  fins,  are  actually  gross  exaggerations  arising 

from  two  possible  phenomena.  The  first  is  the  fact  that  the 

front  fins  are  attacned  to  the  missile  body  only  by  a  thin 

quadrilateral  plate  element,  result ii,g  in  a  lack  of  stiffness 

for  these  fins.  The  second  is  the  fact  thai  scaling  for  each 

mode  shape  in  a  normal  modes  analysis  is  ai'bicrary  since  no 

external  load  is  applied;  only  the  relative  shape  is 

important.  The  deformations  of  the  front  fins  in  those  modes 

was  actually  less  than  0.3  inche.s . 


Figure  17  AlM-9  Missile  Model  Mode  4  Response 


Figure  19  AIM-9  Missile  Model  Mode 


Response 
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APPENDIX  C 


TRADE-OFF  STUDIES 

Several  studies  were  conducted  on  the  basic  AIM-9  model  in 
an  attempt  to  determine  the  effects  of  changing  the  location 
of  the  point  masses  and  the  material  properties  of  the 
structure.  The  results  of  these  trade-off  studies  are  given 
in  the  five  cases  discussed  below  in  the  order  in  which  they 
were  conducted,  with  the  current  AIM-9  model  being  Case  V. 
For  each  of  the  cases,  the  first  six  mode  shapes  and  their 
corresponding  natural  frequencies  are  described,  except  for 
Case  II,  in  which  only  the  first  four  mode  shapes  were 
determined.  In  studying  these  cases,  one  can  see  that 
changing  the  mass  location  or  material  properties  can  change 
both  the  mode  shapes  and  the  resonant  frequencies,  sometimes 
quite  drastically. 


Case  I  —  missile  body  only  with  point  masses  located  on  top  of 
missile  body. 

•  Mode  1  <50.02  hz)  —  bending  in  x-z  plane 

•  Mode  2  <53.36  hz)  —  bending  in  y-z  plane 

•  Mode  3  <71.39  hz)  —  rotation  about  z— axis 

•  Mode  4  <209.60  hz)  -  complex  bending  in  y-z  plane 

•  Mode  5  <262,70  hz)  -  complex  bending  in  y-z  plane  and  body 
torsion  in  x~y  plane 
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Mode  6  (291.80  hz)  -  complex  bending  in  y-z  plane  and  body, 
torsion  in  x-y  plane 


Case  II  -  missile  body  only  with  masses  located  on  bottom  of 
missile  body  but  in  same  z-axis  location. 

•  Mode  1  (33.42  hz)  -  rotation  about  z-axis,  with  bending  in 
x-z  plane 

•  Mode  2  (53.02  hz)  -  bending  in  y-z  plane 

•  Mode  3  (59.75  hz)  -  bending  in  x-z  plane 

•  Mode  4  (140.60  hz)  -  complex  bending  in  y-z  plane 

Case  III  -  full  AIM-9  model  with  fins  and  wings,  and  five 
point  masses  located  on  bottom  of  missile  body. 

»  Mode  1  (24.56  hz)  -  rotation  about  z-axis  (about  the 

attaclinient  points) 

•  Mode  2  (50  80  hz)  -  bending  in  y-z  plane 

•  Mode  3  (50.85  hz)  —  bending  in  x-z  plane  with  -some  body 

torsion  in  x-y  plane 

•  Mode  4  (57.36  hz)  -  front  fin  rotation  about  z-axis;  no 

body  motion 

•  Mode  5  (57.72  hz)  -  front  fin  rotation  about  z-axis  with 
some  body  torsion  in  x-y  plane 

•  Mode  6  (58.70  hz)  -  front  fin  rotation  about  z-axis  with 

some  body  torsion  in  x-y  plane  and  bending  in  x-z  plane 

Case  IV  -  full  AIM-9  model  with  ten  point  masses  located  on 
top  and  bottom  of  m.issile  body  at  the  same  z-axis  location. 

•  Mode  1  (26,68  hz)  -  rotation  about  z-axis  (about  the 

attachment  points) 

•  Mode  2  (50.58  hz)  -  bending  in  x-z  plane 

•  Mode  3  (51.01  hz)  -  bending  in  y-z  plane 
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•  Mode  4  (57.36  hz)  -  front  fin  rotation  about  z-axis;  no 

body  motion 

•  Mode  5  (57.57  hz)  -  front  fin  rotation  about  z-axis  with 

some  body  torsion  in  x-y  plane 

•  Mode  6  (58.71  hz)  -  front  fin  rotation  about  z-axis  with 
some  body  torsion  in  x-y  plane  and  bending  in  x-z  plane 

Case  V  -  full  AIM-9  model  with  five  point  masses  located  on 
bottom  of  missile  body  and  the  rear  wings  changed  from  steel 
to  aluminum. 

•  Modes  given  previously  in  Appendix  B. 
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APPENDIX  D 

MODAL  TEST  SYSTEM  TEST  EQUIPMENT 

Scientific-Atlanta  Spectrviin  Analyzer  Model  SD-380 

Tne  Model  SD-380  is  a  microprocessor-based  one,  two  or 
four  channel  signal  analyzer.  It  measures  the  signal (s) 
present  at  one,  two  or  four  BNC  connectors  located  on  the  rear 
panel,  and  in  the  case  of  the  two  and  four  channel  options, 
the  relationship  betvreen  two  signals.  In  addition  to  standard 
signal  analysis  capabilities,  it  is  equipped  with  a  built-in 
waoerfall  (or  cascade)  feature,  a  56K  word  input  memory,  an 
IEEE  interface,  and  the  ability  to  perform  single/double 
dif terentiation/integration .  Adaitional  specifications  and 
information  on  the  SL)-380  are  supplied  by  Scientif ic-Atlanta 
[Ref.  14]  . 

COMPAQ  Deskpro  386/25e  Model  60  Computer  with  EMODAL  Software 

The  Deskpro  386/25e  includes  a  25  Mhz  386  microprocessor, 
four  (Mb)  of  random  access  memory,  a  60  Mb  hard  disk  drive,  a 
high-density  1.2  Mb,  5  1/4  inch  floppy  di.sk  drive,  and  a  high- 
density  1.44  Mb,  3  1/2  inch  floppy  disk  drive.  It  comes 

equipped  with  an  enhanced  keyboard,  a  COMPAQ  Video  Graphics 
monitor,  and  a  mouse. 
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The  software  used  in  this  test  setup  is  Entek  Scientific 
Corporation's  EMODAL  modal  analysis  software  package.  It  can 
he  used  as  a  stand-alone  modal  testing  tool,  or  as  a  tool  to 
collect  and  process  data.  It  has  been  designed  to  run  with  or 
without  peripheral  devices;  however,  a  spectrum  analyzer  must 
be  connected  when  acquiring  measurement  data.  To  allow  the 
EMODAL  software  to  communicate  with  the  spectrum  analyzer,  tlie 
Deskpro  386/25e  has  had  a  National  Instruments  AT-GPIB  board 
installed,  which  supplies  tht  necessary  IEEE-488  interface. 
More  specific  information  on  the  EMODAL  software  can  be 
obtained  from  the  EMODAL-PC  Instructioii  Manual  [Ref.  3]  . 

Bruel  &  Kjaer  (B&K)  Exciter  Tyi-e  4801 

The  Type  4801  exciter  was  designed  to  be  combined  with  anv 
one  of  several  different  exciter  heads,  resulting  in  a 
complete  exciter  (or  shaker)  assembly,  capable  of  generating 
up  to  a  100  pound  excitation  force.  It  is  not  considered  a 
portable  shaker  with  its  rather  substantial  weight  of  180 
pounds,  but  it  can  be  positioned  at  any  angle  to  provide 
optimum  flexibility  and  use.  With  the  general  purpose  head, 
the  Type  4801  exciter  is  capable  of  providing  displacements  up 
to  0.5  inches.  It  has  an  internal  natural  frequency  of  7200 
hz,  with  the  base  having  a  natural  frequency  of  10-14  hz . 
More  detailed  information  can  be  found  in  tne  manufacturer's 
Instructions  and  Applications  booklet  [Ref.  IG]  and  Vibration 
Exciter  System  V  Instruction  Manual  [Ref.  17] . 
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B&K  General  Purpose  £xcit.er  Head  Type  4812 

The  Type  4812  general  purpose  head  was  designed  to  provide 
a  low  distortion  operating  band  from  DC  to  10,000  hz  when  used 
with  the  Type  4801  exciter  and  the  Type  2707  amplifier.  The 
'g'  loads  available  for  testing  vary  with  the  test  frequency 
and  run  from  a  low  of  zero  in  the  DC  range  to  lOOg  at  150-160 
hz .  This  general  purpose  impedance  head  has  both  a  force 
transducer  and  a  response  transducer  incorporated  into  it,  and 
is  used  to  obtain  driving  point  measurements .  More  detailed 
information  can  be  found  in  the  manufacturer's  Instructions 
and  Applications  booklet  [Ref.  18] . 

B&K  Power  Amplifier  Type  2707 

The  Type  2707  power  amplifier  was  designed  to  provide 
proper  power  and  protection  to  the  Type  4801  exciter. 
Protective  circuits  incorporated  into  the  power  amplifier 
include  power  phase  protection,  signal  ground  fault, 
overtemperature,  displacement,  overcurrent  and  waveform 
clipping  indications.  More  detailed  information  can  be  found 
in  the  manufacturer's  Instruction  Manual  for  Power  Amplifier 
Type  2707  [Ref.  19] . 

PCB  Piazoelectronics  Flexcel  Modal  Accelerometor  Model  336A 

The  Model  336A  modal  accelerometer  functions  to  transfer 
vibratory  motion  into  a  measurable  voltage  signal  for 
convenient  display  and  analysis.  It  is  a  low-mass 


accQleromcter  which  s«3rv'a.s  to  minimico  mass-loading  effects. 
Cliaracteri.stics  of  this  acceleiometei  important  in  modal 
analysis  applications  include  high  output,  low  noise,  matched 
phase  response,  low  transient  thermal  sensitivity,  and  most 
importantly  for  the  missile  modal  test,  good  low  frequency 
response,  provided  by  a  composite  construction  which  isolates 
from  ground  and  thermal  transients.  Since  this  accelerometer 
contains  a  built-in  microelectronic  charge-to-voltage 
amplifier,  a  constant  currei't  power  unit  is  required  to 
of'erate  it  .  Specifics  concerning  the  Model  336A  include  a 
frequency  raiigu  of  1-2000  hi:,  a  sensitivity  of  1000  mV/g,  and 
a  resor.aiit  frequency  (mounted)  of  greater  than  9  khz .  More 
detailed  information  can  be  found  .in  the  manufacturer's 
operating  instructions  [Ref.  20]. 

PCB  Piezoelectronics  Dual-Mocie  Charge  Amplifier  Model  464A 

The  Model  4  64A  is  a  dxial-mode  device  which  can  be  operated 
either  as  a  charge  amplifier  (for  transducers  without  builu-in 
microelectronics)  or  as  an  Integrated-Circuit-Piezoelectric 
(ICP)  constant  current  power  source  (for  transducers  with 
built-in  microelectronics) .  Operation  is  simple  in  that 
switching  to  the  mode  of  operation  (charge  or  ICP) ,  dialing  in 
the  transducer  sensitivity,  and  switching  to  the  desired  range 
is  all  that  is  required.  The  Model  464A  has  a  frequency 
response  from  near  DC  to  100  khz.  For  this  test  setup,  a 
plug-in  180  khz  low-pas.s  filter  was  installed,  as  this  was  all 
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that  was  available  at  the  time  oi  the  test  setup.  howev«i.- 
the  tilters  are  available  ii-  varying  sizes,  rangiiig  troni  1  khz 
to  180  khz.  More  detailed  information  can  be  found  in  the 
manutr.ctui  er' s  operating  instructions  LRuf  .  21]. 

PCfi  Piezoelactronica  Mechanical  Impadanca  Sansor  Modal  288A11 

The  Model  ISSAll  is  a  low  mass,  combination  force 
transducer  and  accelerometer  which  functions  to  transfer  Loth 
the  excitation  (force)  and  the  response  (motion)  into 
electrical  signals.  The  built-in  microelectronic  amplifier 
converts  the  high  impedance  signals  generated  by  the  quartz 
elemei^t  force  gauge  and  accelerometer  to  relatively  low 
impedance  vortage  signals.  Since  this  sensor  contains  a 
built-in  amplifier,  a  constant,  current  power  unit  is  required 
to  operate  it.  Cliaracteristics  of  this  sensor  include  low 
mass  (only  45  grams) ,  nominal  sensitivity  of  100  mV/g  (motion) 
and  1000  mV/g  (force),  and  frequency  range  of  5-5000  hz .  More 
detailed  information  can  be  found  in  the  manufacturer's 
operating  instructions  [Ref.  22] . 
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